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Abstract. We give a survey of algorithms for computing topological invari- 
ants of semi-algebraic sets with special emphasis on the more recent devel- 
opments in designing algorithms for computing the Betti numbers of semi- 
algebraic sets. Aside from describing these results, we discuss briefly the back- 
ground as well as the importance of these problems, and also describe the main 
tools from algorithmic semi-algebraic geometry, as well as algebraic topology, 
which make these advances possible. We end with a list of open problems. 
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1. Introduction 





This article has several goals. The primary goal is to provide the reader with 
a thorough survey of the current state of knowledge on efficient algorithms for 
computing topological invariants of semi-algebraic sets - and in particular their 
Betti numbers. At the same time we want to provide graduate students who intend 
to pursue research in the area of algorithmic semi-algebraic geometry, a primer on 
the main technical tools used in the recent developments in this area, so that they 
can start to use these themselves in their own work. Lastly, for experts in closely 
related areas who might want to use the results described in the paper, we want to 
present self-contained descriptions of these results in a usable form. 

With this in mind we first give a short introduction to the main algorithmic 
problems in semi-algebraic geometry, their history, as well as brief descriptions of 
the main mathematical and algorithmic tools used in the design of efficient algo- 
rithms for solving these problems. We then provide a more detailed description of 
the more recent advances in the area of designing efficient algorithms for computing 
the Betti numbers of semi-algebraic sets. Since the design of these algorithms draw 
on several new ideas from diverse areas, we describe some of the most important 
ones in some detail for the reader's benefit. The goal is to provide the reader with a 
short but comprehensive introduction to the mathematical tools that have proved 
to be useful in the area. The reader who is interested in being up-to-date with 
the recent developments in this area, but not interested in pursuing research in the 
area, can safely skip the more technical sections. Throughout the survey we omit 
most proofs referring the reader to the appropriate references where such proofs 
appear. 

The rest of the paper is organized as follows. In Section |2] we discuss the back- 
ground, significance, and history of algorithmic problems in semi-algebraic geome- 
try and topology. In Section [3] we state some of the recent results in the field. In 
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Section |4] we outline a few of the basic algorithmic tools used in the design of algo- 
rithms for dealing with semi-algebraic sets. These include the cylindrical algebraic 
decomposition, as well as the critical point method exemplified by the roadmap 
algorithm. In Section [5] we provide the reader some relevant facts and definitions 
from algebraic topology which are used in the more modern algorithms, including 
definitions of cohomology of simplicial complexes as well as semi-algebraic sets, the 
Nerve Lemma and its generalizations for non-Leray covers, the descent spectral 
sequence and the basic properties of homotopy colimits. In Section [6] we describe 
recent progress in the design of algorithms for computing the higher Betti numbers 
of semi-algebraic sets. In Section [7] we restrict our attention to sets defined by qua- 
dratic inequalities, and describe recent progress in the design of efficient algorithms 
for computing the Betti numbers of such sets. In Section|8]we describe a simplified 
version of an older algorithm for efficiently computing the Betti numbers of an ar- 
rangement - where the emphasis is on obtaining tight bounds on the combinatorial 
complexity only (the algebraic part of the complexity being assumed to be bounded 
by a constant). We end by listing some open problems in Section [9] 
Prerequisites. In this survey we are aiming at a wide audience. We expect that 
the reader has a basic background in algebra, has some familiarity with simplicial 
complexes and their homology, and the theory of NP and ^/^P-completeness. Beyond 
these we make no additional assumption of any prior advanced knowledge of semi- 
algebraic geometry, algebraic topology, or the theory of computational complexity. 



2. Semi-algebraic Geometry: Background 

2.1. Notation. We first fix some notation. Let R be a real closed field (for exam- 
ple, the field R of real numbers or Maig of real algebraic numbers). A semi-algebraic 
subset of R'' is a set defined by a finite system of polynomial equalities and inequali- 
ties, or more generally by a Boolean formula whose atoms are polynomial equalities 
and inequalities. Given a finite set V of polynomials in R[Xi, . . . , X^], a subset S 
of R*^ is T'-semi-algebraic if S is the realization of a Boolean formula with atoms 
P = 0,P>OorP<0 with P G V. It is clear that for every semi-algebraic subset 
S of R'^ there exists a finite set V of polynomials in R[Xi, . . . , Xk] such that S is 
T'-semi-algebraic. We call a semi-algebraic set a P-closed semi-algebraic set if it 
is defined by a Boolean formula with no negations with atoms P = 0, P > 0, or 
P < with P eV. 

For an element a e R we let 

{0 if a = 0, 
1 if a > 0, 
-1 ifa<0. 

A sign condition on V is an element of {0,1,-1}^. For any semi-algebraic set 
Z C R'" the realization of the sign condition a over Z, TZ{<7, Z), is the semi-algebraic 
set 

{xeZ \ /\ sign(P(x)) = a(P)}, 
Pev 

and in case Z = R*^ we will denote 7?.(ct, Z) by just TZ{a). 
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If 7^ is a finite subset of R[Xi, . . . , X^], we write the set of zeros of V in R*^ as 

Z(7',R'=) = {xeR'' \ /\ P{x) = 0}. 

Pev 

We will denote by Bk{0,r) the open ball with center and radius r in R''. We 
will also denote by S*" the unit sphere in R*^^^ centered at the origin. Notice that 
these sets are semi-algebraic. 

For any semi-algebraic set X, we denote by X the closure of X, which is also a 
semi-algebraic set by the Tarksi-Seidenberg principle [501 [53] (see [H] for a modern 
treatment) . The Tarksi-Seidenberg principle states that the class of semi-algebraic 
sets is closed under linear projections or equivalently that the first order theory 
of the reals admits quantifier elimination. It is an easy exercise to verify that 
the closure of a semi-algebraic set admits a description by a quantified first order 
formula. 

For any semi-algebraic set S, we will denote by bi{S) its i-th Betti number, 
which is the dimension of the i-th cohomology group, H*(5', Q), taken with rational 
coefficients, which in our setting is also isomorphic to the i-th homology group. 



Hi(5', Q) (see Section 5.3 below for precise definitions of these groups). In particular, 
bo{S) is the number of semi-algebraically connected components of S. We will 
sometimes refer to the sum b{S) = bi{S) as the topological complexity of a semi- 

i>0 

algebraic set S. 

Remark 2.1. Departing from usual practice, in the description of the algorithms 
occurring later in this paper we will mostly refer to the cohomology groups instead 
of the homology groups. Even though the geometric interpretation of the coho- 
mology groups is a bit more obscure than that for homology groups (see Section 
5.1.1 below), it turns out that from the point of view of designing algorithms for 
computing Betti numbers of semi-algebraic sets (at least for those discussed in this 
survey) the usual geometric interpretation of homology as measuring the number 
of "holes" or "tunnels" etc. is of little use, and the main concepts behind these 
algorithms are better understood from the cohomological point of view. This is the 
reason why we emphasize cohomology over homology in what follows. 

2.2. Main Algorithmic Problems. Algorithmic problems in semi-algebraic ge- 
ometry typically consist of the following. We are given as input a finite family, 
V C R[Ai, . . . ,Xk], as well as a formula defining a 7^-semi-algebraic set S. The 
task is to decide whether certain geometric and topological properties hold for S, 
and in some cases also computing certain topological invariants of S. Some of the 
most basic problems include the following. 
Given a P-semi- algebraic set C R*": 

(1) decide whether it is empty or not, 

(2) given two points x,y G S, decide if they are in the same connected compo- 
nent of S and if so output a semi-algebraic path in S joining them, 

(3) compute semi-algebraic descriptions of the connected components of S, 

(4) compute semi-algebraic descriptions of the projection of S onto some linear 
subspace of R*^ (this problem is also known as the quantifier elimination 
problem for the first order theory of the reals and many other problems can 
be posed as special cases of this very general problem) . 
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At a deeper level we have problems of more topological flavor, such as: 
(5) compute the cohomology groups of S, its Betti numbers, its Euler-Poincare 
characteristic etc.. 



(6) compute a semi-algebraic triangulation of 5* (cf. Definition 4.4 below), as 
well as 

(7) compute a decomposition of S into semi-algebraic smooth pieces of various 
dimensions which fit together nicely (a Whitney-regular stratification) . 

The complexity of an algorithm for solving any of the above problems is measured 
in terms of the following three parameters: 

• the number of polynomials, s = =ffV, 

• the maximum degree, d — maxpg-p deg(P), and 

• the number of variables, k. 

Definition 2.2 (Complexity). A typical input to the algorithms considered in this 
survey will be a set of polynomials with coefficients in an ordered ring D (which 
can be taken to be the ring generated by the coeffcients of the input polynomials). 
By complexity of an algorithm we will mean the number of arithmetic operations 
(including comparisons) performed by the algorithm in the ring D. In case the 
input polynomials have integer coefficients with bounded bit-size, then we will 
often give the bit-complexity, which is the number of bit operations performed by 
the algorithm. We refer the reader to L22j [Chapter 8] for a full discussion about the 
various measures of complexity. 

Even though the goal is always to design algorithms with the best possible 
complexity in terms of all the parameters s, d, k, the relative importance of the 
parameters is very much application dependent. For instance, in applications in 
computational geometry it is the combinatorial complexity (that is the dependence 
on s) that is of paramount importance, the algebraic part depending on d, as well 
as the dimension k, are assumed to be bounded by constants. On the other hand in 
algorithmic real algebraic geometry, and in applications in complexity theory, the 
algebraic part depending on d is considered to be equally important. 

2.3. Brief History. Even though there exist algorithms for solving all the above 
problems, the main research problem is to design efficient algorithms for solving 
them. The complexity of the first decision procedure given by Tarski |60] to solve 
Problems 1 and 4 listed in Section |2.2| is not elementary recursive, which implies 
that the running time cannot be bounded by a function of the size of the input 
which is a fixed tower of exponents. The first algorithm with a significantly better 
worst-case time bound was given by Collins ^34] in 1976. His algorithm had a worst 
case running time doubly exponential in the number of variables. Collins' method 
is to obtain a cylindrical algebraic decomposition of the given semi-algebraic set 



(see Section 4.1 below for definition). Once this decomposition is computed most 
topological questions about semi-algebraic sets such as those listed in Section 2.2 can 
be answered. However, this method involves cascading projections which involve 
squaring of the degrees at each step resulting in a complexity which is doubly 
exponential in the number of variables. 

Most of the recent work in algorithmic semi-algebraic geometry has focused on 
obtaining single exponential time algorithms - that is algorithms with complexity 
of the order of (sc?)'^°'^' rather than (sd)'^'' . An important motivating reason be- 
hind the search for such algorithms, is the following theorem due to Gabrielov and 
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Vorobjov jlQ] (see [56l |65l ESI E] , as well as the survey article [21], for work leading 
up to this result). 

Theorem 2.3. [JD] For a V -semi-algebraic set S C R*^, the sum of the Betti 
numbers of S (refer to Section^ below for definition) is bounded by {0{s'^d))^ , 
where s — and d — maxpgp deg(P). 

For the special case of 7^-closed semi-algebraic sets the following slightly better 
bound was known before (and this bound is used in an essential way in the proof 



of Theorem 2.3 1. Using the same notation as in Theorem 2.3 above we have 



Theorem 2.4. [5J For a V-closed semi- algebraic set S C R*^, the sum of the Betti 
numbers of S is bounded by (O(sd))'^. 

Remark 2.5. These bounds are asymptotically tight, as can be already seen from 
the example where each P £ V is a, product of d generic polynomials of degree one. 
The number of connected components of the T'-semi-algebraic set defined as the 
subset of R*^ where all polynomials in V are non-zero is clearly bounded from below 
by (nisd))". 

Notice also that the above bound has single exponential rather than double expo- 
nential dependence on k. Algorithms with single exponential complexity have now 
been given for several of the problems listed above and there have been a sequence 
of improvements in the complexities of such algorithms. We now have single ex- 
ponential algorithms for deciding emptiness of semi-algebraic sets [42) , 143 ] l57 l [T6]. 
quantifier elimination [STJ [111 M i deciding connectivity [SHI Sll 133 ISHl [IZ] , com- 
puting descriptions of the connected components [371 H^i computing the Euler- 



Poincare characteristic (see Section 5.3.1 below for definition) [51 jTH], as well as 



the first few (that is, any constant number of) Betti numbers of semi-algebraic sets 
[20l [To] . These algorithms answer questions about the semi-algebraic set S with- 



out obtaining a full cylindrical algebraic decomposition (see Section 4.1 below for 
definition) , which makes it possible to avoid having double exponential complexity. 
Moreover, polynomial time algorithms are now known for computing some of these 
invariants for special classes of semi-algebraic sets [2 HSl [3 El IH] . We describe 
some of these new results in greater detail in Section [3] 

2.4. Certain Restricted Classes of Semi-algebraic Sets. Since general semi- 



algebraic sets can have exponential topological complexity (cf. Remark 2.5), it is 
natural to consider certain restricted classes of semi-algebraic sets. One natural 
class consists of semi-algebraic sets defined by a conjunction of quadratic inequali- 
ties. 

2.4.1. Quantitative Bounds for Sets Defined by Quadratic Inequalities. Since sets 
defined by linear inequalities have no interesting topology, sets defined by quadratic 
inequalities can be considered to be the simplest class of semi-algebraic sets which 
can have non-trivial topology. Such sets are in fact quite general, since every semi- 
algebraic set can be defined by a (quantified) formula involving only quadratic 
polynomials (at the cost of increasing the number of variables and the size of the 
formula). Moreover, as in the case of general semi-algebraic sets, the Betti numbers 
of such sets can be exponentially large. For example, the set S" C R'^ defined by 

X,{l-X,)<0,...,Xk{l-Xk)<0, 
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has boiS) ^ 2^ 

Hence, it is somewhat surprising that for any constant ^ > 0, the Betti num- 
bers bk-i{S), . . . , bk-i{S), of a basic closed semi-algebraic set S* C R'^ defined by 
quadratic inequalities, are polynomially bounded. The following theorem which ap- 
pears in [7] is derived using a bound proved by Barvinok U on the Betti numbers 
of sets defined by few quadratic equations. 

Theorem 2.6. [7 Let R a real closed field and S C R'^ be defined by 
Pi<0,...,Ps<0, deg(P,) < 2, 1 < i < s. 

Then, for any £ > 0, 

Notice that for fixed £ this gives a polynomial bound on the highest £ Betti 
numbers of S (which could possibly be non-zero). Observe also that similar bounds 
do not hold for sets defined by polynomials of degree greater than two. For instance, 
the set Van'' defined by the single quartic inequality, 

k 

J2xhX^-lf-e>0, 

i=l 

will have bk^i{V) — 2*^, for all small enough e > 0. 

To see this observe that for all sufficiently small e > 0, R*" \ is defined by 

A; 

4=1 

and has 2'^ connected components since it retracts onto the set {O,!}*^. It now 
follows that 

&fe-l("^^) = foo(R'\^^) = 2^ 

where the first equality is a consequence of the well-known Alexander duality the- 
orem (see [Ml pp. 296]). 

2.4.2. Relevance to Computational Complexity Theory. Semi-algebraic sets defined 
by a system of quadratic inequalities have a special significance in the theory of com- 
putational complexity. Even though such sets might seem to be the next simplest 
class of semi-algebraic sets after sets defined by linear inequalities, from the point 
of view of computational complexity they represent a quantum leap. Whereas there 
exist (weakly) polynomial time algorithms for solving linear programming, solving 
quadratic feasibility problem is provably hard. For instance, it follows from an easy 
reduction from the problem of testing feasibility of a real quartic equation in many 
variables, that the problem of testing whether a system of quadratic inequalities is 
feasible is NPpt-complete in the Blum-Shub-Smale model of computation (see j29j ) . 
Assuming the input polynomials to have integer coefficients, the same problem is 
NP-hard in the classical Turing machine model, since it is also not difficult to see 
that the Boolean satisfiability problem can be posed as the problem of deciding 
whether a certain semi-algebraic set defined by quadratic inequalities is empty or 
not. 

Counting the number of connected components of such sets is even harder. In 
fact, it is shown in ,11] that for £ < log/c, computing the ^-th Betti number of 
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a basic semi-algebraic set defined by quadratic inequalities in R'^ is ^P-hard. In 
contrast to these hardness results, the polynomial bound on the top Betti numbers 
of sets defined by quadratic inequalities gives rise to the possibility that these might 
in fact be computable in polynomial time. 

2.4.3. Projections of Sets Defined by Few Quadratic Inequalities. A case of inter- 
mediate complexity between semi-algebraic sets defined by polynomials of higher 
degrees and sets defined by a fixed number of quadratic inequalities is obtained by 
considering linear projections of such sets. The operation of linear projection of 
semi-algebraic sets plays a very significant role in algorithmic semi- algebraic geom- 
etry. It is a consequence of the Tarski-Seidenberg principle (see for instance [HJ 
Theorem 2.80]) that the image of a semi-algebraic set under a linear projection is 
semi-algebraic, and designing efficient algorithms for computing properties of pro- 
jections of semi-algebraic sets (such as its description by a quantifier-free formula) 
is a central problem of the area and is a very well-studied topic (see for example 
[571 [ini E] or [551 Chapter 14]). However, the complexities of the best algorithms 
for computing descriptions of projections of general semi-algebraic sets is single ex- 
ponential in the dimension and do not significantly improve when restricted to the 
class of semi-algebraic sets defined by a constant number of quadratic inequalities. 
Indeed, any semi-algebraic set can be realized as the projection of a set defined by 
quadratic inequalities, and it is not known whether quantifier elimination can be 
performed efficiently when the number of quadratic inequalities is kept constant. 
However, it is shown in [24j that, with a fixed number of inequalities, the projections 
of such sets are topologically simpler than projections of general semi-algebraic sets. 

More precisely, let S C M*^^™ be a closed and bounded semi-algebraic set defined 
by Pi > 0, . . . , > 0, with P, e . . . , Xfc, Fi, . . . , r™],deg(P,) < 2, 1 < z < £ 

(For technical reasons, which we do not delve into, it is necessary in this case to 
restrict ourselves to the case where R = K.) Let tt : M*''^™ —i- be the projection 
onto the last m coordinates. In what follows, the number of inequalities, £, used 
in the definition of S will be considered as fixed. Since, tt{S) is not necessarily 
describable using only quadratic inequalities, the bound in Theorem |2.6| does not 
hold for it{S) and n{S) can in principle be quite comphcated. Using the best known 
complexity estimates for quantifier elimination algorithms over the reals (see [551 
Chapter 14]), one gets single exponential (in k and m) bounds on the degrees and 
the number of polynomials necessary to obtain a semi-algebraic description of tt{S). 
In fact, there is no known algorithm for computing a semi-algebraic description of 
it{S) in time polynomial in k and m. Nevertheless, we know that for any constant 
q > 0, the sum of the first q Betti numbers of tt{S) is bounded by a polynomial in 
k and m. 

Theorem 2.7. [51j Let S C R'^^™' be a closed and bounded semi-algebraic set 
defined by 

Pi > 0, . . . , P, > 0, P, e M[Xi, . . . , Xfc, n, . . . , y,„], deg(P,) < 2, 1 < i < £ 

Let TT : M'^^™ —f M™ be the projection onto the last m coordinates. For any q > 
0, 0<q<k, 

(2.1) J2h{n{S))<{k + mf^'^'\ 

i=0 
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This suggests, from the point of view of designing efficient (polynomial time) 
algorithms in semi- algebraic geometry, that images under linear projections of semi- 
algebraic sets defined by a constant number of quadratic inequalities, are simpler 
than general semi-algebraic sets. So they should be the next natural class of sets 
to consider, after sets defined by linear and quadratic inequalities. 

2.5. Some Remarks About the Cohomology Groups. Since in this survey 
we focus mainly on the algorithmic problem of computing the Betti numbers of 
semi-algebraic sets, which are the dimensions of the various cohomology (also ho- 
mology) groups of such sets, it is perhaps worthwhile to say a few words about our 
motivations behind computing them, and also their connections with other parts 
of mathematics, especially with computational complexity theory. 

2.5.1. Motivation behind computing the zero-th Betti number. The algorithmic prob- 
lems of deciding whether a given semi-algebraic set is empty or if it is connected, 
have obvious applications in many different areas of science and engineering. (Re- 
call that the number of connected components of a semi-algebraic set S is equal to 
its zero-th Betti number, hQ{S).) For instance, in robotics, the configuration space 
of a robot can be modeled as a semi-algebraic set. Similarly, in molecular chemistry 
the conformation space of a molecule with constraints on bond lengths and angles 
is a semi-algebraic set. In both these cases understanding connectivity information 
is important: for solving motion planning problem in robotics, or for determining 
possible molecular conformations in molecular chemistry. 

2.5.2. The higher Betti numbers. The higher cohomology groups of semi-algebraic 
sets, which measure higher dimensional connectivity, do not appear to have such 
obvious applications. Nevertheless, there exist several reasons why the problem 
of computing the higher homology groups of semi-algebraic sets is an important 
problem and we mention a few of these below. 

Firstly, the algorithmic problem of pinning down the exact topology of any given 
topological space, such as a semi-algebraic set in M*^, is an exceedingly difficult 
problem. In fact, the general problem of determining if two given spaces are home- 
omorphic is undecidable |5Q] . In order to get around this difficulty, mathematicians 
since the time of Poincare have devised more easily computable (albeit weaker) 
invariants of topological spaces. One reason that cohomology groups are so im- 
portant, is that unlike other topological invariants, they are readily computable - 
they allow one to discard a large amount of information regarding the topology of a 
given space, while retaining just enough to derive important qualitative information 
about the space in question. For instance, in the case of semi-algebraic sets, the 
dimensions of the cohomology groups also known as the Betti numbers, determine 
qualitative information about the set, such as connectivity (in the usual sense), 
number of holes and/or tunnels (i.e. higher dimensional connectivity), its Euler- 
Poincare characteristic (a discrete valuation with properties analogous to those of 
volume) etc. 

Secondly, the reach of cohomology theory is not restricted to the continuous 
domain (such as the study of algebraic varieties in C*^ or semi-algebraic sets in M'^). 
As a consequence of a series of astonishing theorems (conjectured by Andre Weil 
[68] and proved by Deligne [36l[37], Dwork [38] et al.), it turns out that the number 
of solutions of systems of polynomial equations over a finite field, F^, in algebraic 
extensions of F^, is governed by the dimensions of certain (appropriately defined) 
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cohomology groups of the associated variety (see below) . In this way, cohomology 
theory plays analogous roles in the discrete and continuous settings. 

Finally, the algorithmic problem of computing the cohomology groups of semi- 
algebraic sets is important from the viewpoint of computational complexity theory 
because of the following. It is easily seen that the classical NP-complete problem 
in discrete complexity theory, the Boolean satisfiability problem, is polynomial 
time equivalent to the problem of deciding whether a given system of polynomial 
equations in many variables over a finite field (say Z/2Z) has a solution. The real (as 
well as the complex) analogue of this problem has been proved to be NP-complete 
in the real (resp. complex) version of Turing machines, namely the Blum-Shub- 
Smale machine (see [29j). The algebraic variety defined by a system of polynomial 
equations clearly has further structure apart from being merely empty or non-empty 
as a set. In the discrete case, we might want to count the number of solutions - 
and this turns out to be a #P-complete problem. Recently, a #P-completeness 
theory has been proposed for the BSS model as well f271 HSj ~ and the natural #P 
complete problem in this context is computing the Euler-Poincare characteristic of 
a given variety (the Euler-Poincare characteristic being a discrete valuation is the 
"right" notion of cardinality for infinite sets in this context). 

If one is interested in more information about the variety, then in the discrete case 
one could ask to count the number of solutions of the given system of polynomials 
not just over the ground field Fg, but in every algebraic extension, Fgn of the 
ground field. Even though this appears to be an infinite sequence of numbers, its 
exponential generating function (the so called zeta-function of the variety) turns 
out to be a rational function (conjectured by Weil [55], and proved by Dwork [55] ) 
of the form, 

^ Plit)P3{t)...P,rn-l{t) 
^ ' P2{t)P4{t)...P2m{t) ' 

where each Pi is a polynomial with coefficients in a field of characteristic 0, and the 
degrees of the polynomials Pi{t) are the dimensions of (appropriately defined) coho- 
mology groups associated to the variety V defined by the given system of equations. 
In the real and complex setting, the ordinary topological Betti numbers are con- 
sidered some of the most important computable invariants of varieties and carry 
important topological information. Thus, the algorithmic problem of computing 
Betti numbers of constructible sets or varieties, is a natural extension of some of 
the basic problems appearing in computational complexity theory - namely de- 
ciding whether a given system of polynomial equation is satisfiable, and counting 
the number of solutions. This is true in both the discrete and continuous settings. 
Even though, in this survey we concentrate on the latter, some of the techniques 
developed in this context conceivably have applications in the discrete case as well. 

Also note that, by considering a complex variety C C*" as a real semi-algebraic 
set in K^'^, all results discussed in this survey extend directly (with the same asymp- 
totic complexity bounds) to the corresponding problems (of computing the Betti 
numbers) for complex algebraic varieties, and more generally for constructible sub- 
sets of C'^'. 

3. Recent Algorithmic Results 

In this section we list some of the recent progress on the algorithmic problem of 
determining the Betti numbers of semi-algebraic sets. 
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In |.20^, an algorithm with single exponential complexity is given for com- 



puting the first Betti number of semi-algebraic sets (see Section 6.2 below). 
Previously, only the zero-th Betti number (i.e. the number of connected 
components) could be computed in single exponential time. Another im- 
portant result contained in this paper is the homotopy equivalence between 
an arbitrary semi-algebraic set, and a closed and bounded one (which is 
defined using infinitesimal perturbations of the polynomials defining the 
original set) obtained by a construction due to Gabrielov and Vorobjov 
130]. It was conjectured in [40 that these sets are homotopy equivalent. 
This result is important by itself since it allows, for instance, a single ex- 
ponential time reduction of the problem of computing Betti numbers of 
arbitrary semi-algebraic sets to the same problem for closed and bounded 
ones. 

The above result is generalized in [10], where a single exponential time 
algorithm is given for computing the first f. Betti numbers of semi-algebraic 



sets, where £ is allowed to be any constant (see Section 6.3.1 below). More 
precisely, an algorithm is described that takes as input a description of a 
"P-semi-algebraic set S* C R*^ , and outputs the first I + 1 Betti numbers 
of 5, &o('S'), • • • ,bi{S). The complexity of the algorithm is (srf)*^"*'^' , where 
s = H^i'P) and d — maxpg-p deg(P), which is single exponential in k for H. 
any constant. 

In a polynomial time algorithm is given for computing a constant num- 
ber of the top Betti numbers of semi-algebraic sets defined by quadratic in- 
equalities. If the number of inequalities is fixed then the algorithm computes 



all the Betti numbers in polynomial time (see Section 7.4 below). More pre- 
cisely, an algorithm is described which takes as input a semi-algebraic set, 
S, defined by Pi > 0, . . . , > 0, where each Pi € R[Xi, . . . , Xk] has degree 
< 2, and computes the top I Betti numbers of 5*, hk-i{S), . . . , bk-e{S), in 
polynomial time. The complexity of the algorithm is X^i^o (^^^/c^^*"""'*'"" . 

For fixed i, the complexity of the algorithm can be expressed as s^+^fc^"'*' , 
which is polynomial in the input parameters s and k. For fixed s, we obtain 
by letting i — k, an algorithm for computing all the Betti numbers of S 
whose complexity is . 

In [S] , an algorithm is described which takes as input a closed semi-algebraic 
set, 5" C R'', defined by 

Pi > 0, . . . ,P, > 0,P, e R[Xi,. . .,Xk],deg{P,) < 2, 



and computes the Euler-Poincare characteristic of S (see Section 7.3 below). 
The complexity of the algorithm is k'-"-^\ Previously, algorithms with the 
same complexity bound were known only for testing emptiness (as well as 
computing sample points) of such sets [3J |35] . 

In [24 , a polynomial time algorithm is obtained for computing a con- 
stant number of the lowest Betti numbers of semi-algebraic sets defined 
as the projection of se mi-a lgebraic sets defined by few by quadratic in- 
equalities (see Section 7.5 below). More precisely, let S C M*^^™ be a 



closed and bounded semi-algebraic set defined by Pi > 0, . . . , P^ > 0, where 
Pj e M[Xi, . . . , Xfc, Yi, . . . , Ym], and deg(P,) < 2, 1 < i < £. Let tt denote 
the standard projection from E'^^™ onto M™. An algorithm is described 
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for computing the the first q Betti numbers of 77(5*), whose complexity is 
(fc + to)^ < ' ) ^ fixed q and ^, the bound is polynomial in fc + m. 
• The complexity estimates for all the algorithms mentioned above included 
both the combinatorial and algebraic parameters. As mentioned in Section 
[2j in applications in computational geometry the algebraic part of the com- 
plexity is treated as a constant. In this context, an interesting question is 
how efficiently can one compute the Betti numbers of an arrangement of 
n closed and bounded semi-algebraic sets, Si,. . . ,Sn C R*^, where each Si 
is described using a constant number of polynomials with degrees bounded 
by a constant. Such arrangements are ubiquitous in computational ge- 
ometry (see [T]). A naive approa ch u sing triangulations would entail a 
complexity of 0(n^ ) (see Theorem 4.5 below). This problem is considered 



in [8j where an algorithm is described for computing ^-th Betti number, 

n 

^^(U '^»)' < ^ < fc — 1, using 0(n^+^) algebraic operations. Additionally, 
1=1 

one has to perform linear algebra on integer matrices of size bounded by 
0(n^+^) (see Section [s] below) . All previous algorithms for computing the 
Betti numbers of arrangements triangulated the whole arrangement giving 
rise to a complex of size 0{ti? ) in the worst case. Thus, the complexity of 
computing the Betti numbers (other than the zero-th one) for these algo- 

n 

rithms was 0{n^ ). This is the first algorithm for computing ([^ Si) that 

i=l 

does not rely on such a global triangulation, and has a graded complexity 
which depends on i. 

We should also mention at least one other approach towards computation 
of Betti numbers (of complex varieties) that we do not describe in detail 
in this survey. Using the theory of local cohomology and D-modules, Oaku 
and Takayama [55] and Walther [66l|67], have given explicit algorithms for 
computing a sub-complex of the algebraic de Rham complex of the com- 
plements of complex affine varieties (quasi-isomorphic to the full complex 
but of much smaller size) from which the Betti numbers of such varieties 
as well as their complements can be computed easily using linear algebra. 
For readers familiar with de Rham cohomology theory for differentiable 
manifolds, the algebraic de Rham complex is an algebraic analogue of the 
usual de Rham complex consisting of vector spaces of differential forms. 
The computational complexities of these procedures are not analyzed very 
precisely in the papers cited above. However, these algorithms use Grobner 
basis computations over non-commutative rings (of differential operators) , 
and as such are unlikely to have complexity better than double exponential 
(see |67l Section 2.4]). Also, these techniques are applicable only over al- 
gebraically closed fields, and not immediately useful in the semi-algebraic 
context which is our main interest in this paper, and as such we do not 
discuss these algorithms any further. 

4. Algorithmic Preliminaries 



In this section we give a brief overview of the basic algorithmic constructions 
from semi-algebraic geometry that play a role in the design of more sophisticated 
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algorithms. These include cylindrical algebraic decomposition (Section 4.11, the 



algebraic sets (Section 4.3 1 



critical point method (Section 4.2 1, and the construction of roadmaps of semi- 



4.1. Cylindrical Algebraic Decomposition. As mentioned earlier one funda- 
mental technique for computing topological invariants of semi-algebraic sets is 
through Cylindrical Algebraic Decomposition. Even though the mathematical ideas 
behind cylindrical algebraic decomposition were known before (see for example 
[49, ) , Collins [34 was the first to apply cylindrical algebraic decomposition in the 
setting of algorithmic semi-algebraic geometry. Schwartz and Sharir |58j realized 
its importance in trying to solve the motion planning problem in robotics, as well 
as computing topological properties of semi-algebraic sets. Variants of the basic 
cylindrical algebraic decomposition have also been used in several papers in com- 
putational geometry. For instance in the paper by Chazelle et al. [32' , a truncated 
version of cylindrical decomposition is described whose combinatorial (though not 
the algebraic) complexity is single exponential. This result has found several ap- 
plications in discrete and computational geometry (see for instance [33'). 

Definition 4.1 (Cylindrical Algebraic Decomposition). A cylindrical algebraic de- 
composition of R'^ is a sequence Si, . . . ,Sk where, for each 1 < i < k, Si is a, finite 
partition of R' into semi-algebraic subsets, called the cells of level i, which satisfy 
the following properties: 

• Each cell S G Si is either a point or an open interval. 

• For every 1 < * < ^ and every S € St, there are finitely many continuous 
semi-algebraic functions 

?S,i < • ■ • < ^s/s '■ — * 

such that the cylinder S* x R c R'^^ is the disjoint union of cells of 
which are: 

— either the graph of one of the functions ^s.j, for j — 1, . . . ,is: 

{{x',Xj+i) e 5" X R I Xj+i = £,s,j{x')} , 

— or a band of the cylinder bounded from below and from above by the 
graphs of the functions £,s.j and ^sj+i, for j — 0, . . . where we 
take ,^5,0 = and Ci,£s+i ~ +oo: 

{{x',Xj+i) e S* X R I is,j{x') < Xj+i < ^s,j+i{x')} ■ 
We note that every cell of a cylindrical algebraic decomposition is semi-algebraical- 



ly homeomorphic to an open i-cube (0, 1)* (by convention, (0, 1)*^ is a point). 

A cylindrical algebraic decomposition adapted to a finite family of semi-algebraic 
sets Ti , . . . , is a cylindrical algebraic decomposition of R'^ such that every Ti is 
a union of cells, (see Figure [TJ. 

Definition 4.2. Given a finite set V C R[Xi, . . . ,Xk], a subset S of R'^ is is V- 
invariant if every polynomial P € V has a constant sign (> 0, < 0, or = 0) on 5'. 
A cylindrical algebraic decomposition of R*^ adapted to T' is a cylindrical algebraic 
decomposition for which each cell C Q Sk is 7^-invariant. It is clear that if S is V- 
semi-algebraic, a cylindrical algebraic decomposition adapted to T' is a cylindrical 
algebraic decomposition adapted to S. 
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Figure 1 . Example of cylindrical algebraic decomposition of 
adapted to a sphere. 

One important result which underlies most algorithmic applications of cylindrical 
algebraic decomposition is the following (see |22l Chapter 11] for an easily accessible 
exposition) . 

Theorem 4.3. For every finite set V of polynomials in R[Xi, . . . , Xk], there is a 
cylindrical decomposition ofH'' adapted to V. Moreover, such a decomposition can 
he computed in time {sd)'^'^ ' where s = and d = maxp^-p deg(P). 

The cylindrical algebraic decomposition obtained in Theorem |4.3| can in fact be 
refined to give a semi-algebraic triangulation of any given semi- algebraic set within 
the same complexity bound. 

Recall that 

Definition 4.4 (Semi-algebraic Triangulation). A semi-algebraic triangulation of 
a semi-algebraic set 5' is a simplicial complex K together with a semi-algebraic 
homeomorphism from jA'! to S. 

The following theorem states that such triangulations can be computed for any 
closed and bounded semi-algebraic set with double exponential complexity. 

Theorem 4.5. Let S C R'^ be a closed and bounded semi- algebraic set, and let 
Si, ... ,Sq be semi-algebraic subsets of S. There exists a simplicial complex K in 
R*^ and a semi- algebraic homeomorphism h : \K\ — > S such that each Sj is the 
union of images by h of open simplices of K . Moreover, the vertices of K can be 
chosen with rational coordinates. 
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Moreover, if S and each Si are V -semi-algebraic sets, then the semi-algebraic 
triangulation {K,h) can be computed in time (sd)^ * where s = ^f^V and d ~ 
maxpgp deg(P). 

4.2. The Critical Point Method. As mentioned earlier, all algorithms using 
cylindrical algebraic decomposition have double exponential complexity. Algo- 
rithms with single exponential complexity for solving problems in semi-algebraic 
geometry are mostly based on the critical point method. This method was pio- 
neered by several researchers including Grigoriev and Vorobjov |43l 144] , Renegar 
[57] . Canny [3D], Heintz, Roy and Solerno |17], Basu, Pollack and Roy [IS] amongst 
others. In simple terms, the critical point method is nothing but a method for 
finding at least one point in every semi-algebraically connected component of an 
algebraic set. It can be shown that for a bounded nonsingular algebraic hyper- 
surface, it is possible to change coordinates so that its projection to the Xi-axis 
has a finite number of non-degenerate critical points. These points provide at least 
one point in every semi-algebraically connected component of the bounded nonsin- 
gular algebraic hyper-surface. Unfortunately this is not very useful in algorithms 
since it provides no method for performing this linear change of variables. Moreover 
when we deal with the case of a general algebraic set, which may be unbounded or 
singular, this method no longer works. 

In order to reduce the general case to the case of bounded nonsingular alge- 
braic sets, we use an important technique in algorithmic semi-algebraic geometry - 
namely, perturbation of a given real algebraic set in R*" using one or more infinites- 
imals. The perturbed variety is then defined over a non-archimedean real closed 
extension of the ground field - namely the field of algebraic Puiseux series in the 
infinitesimal elements with coefficients in R. 

Since the theory behind such extensions might be unfamiliar to some readers, 
we introduce here the necessary algebraic background referring the reader to |22l 
Section 2.6] for full detail and proofs. 

4.2.1. Infinitesimals and the Field of Algebraic Puiseux Series. 

Definition 4.6 (Puiseux series). A Puiseux series in s with coefficients in R is a 
series of the form 



with fc G Z, i G Z, fli G R, g a positive integer. 

It is a straightforward exercise to verify that the field of all Puiseux series in e 
with coefficients in R is an ordered field. The order extends the order of R, and e is 
an infinitcsimally small and positive, i.e. is positive and smaller than any positive 



Notation 1. The field of Pusisex series in e with coefficients in R contains as a 
subfield, the field of Puiseux series which are algebraic over R[e]. We denote by 
R(£) the field of algebraic Puiseux series in ( with coefficients in R. 

The following theorem is classical (see for example [22j Section 2.6] for a proof). 

Theorem 4.7. The field R(e) is real closed. 



(4.1) 




i>k 



r G R. 
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Definition 4.8 (The lim^ map). When a e R(£) is bounded by an element of R, 
Hm£(a) is the constant term of a, obtained by substituting for e in a. 

Example 4.9. A typical example of the application of the lim map can be seen in 
Figures [2] and [s] below. The first picture depicts the algebraic set Z(Q,R'^), while 
the second depicts the algebraic set Z{Q,'R{0^) (where we substituted a very small 
positive n umb er for C, in o rder to able display this set), where Q and Q are defined 



by Eqn. (4.4| and Eqn. (4.3 1 resp. The algebraic sets Z(g,R'^) and Z(Q,R(C) ) 



are related by 



Z(Q,R3)^lim Z(Q,R(C)'). 



Since we will often consider the semi-algebraic sets defined by the same formula, 
but over different real closed extensions of the ground field, the following notation 
is useful. 

Notation 2. Let R' be a real closed field containing R. Given a semi-algebraic set 
S in R*^, the extension of S to R', denoted Ext(5', R'), is the semi-algebraic subset 
of R''' defined by the same quantifier free formula that defines S. 

The set Ext(S', R') is well defined (i.e. it only depends on the set S and not on 
the quantifier free formula chosen to describe it). This is an easy consequence of 
the transfer principle. 

We now return to the discussion of the critical point method. In order for 
the critical point method to work for all algebraic sets, we associate to a possibly 
unbounded algebraic set Z C a, bounded algebraic set Z' C R(e)'^^^, whose 
semi-algebraically connected components are closely related to those of Z. 

Let Z = Z((5,R'') and consider 

Z' = Z(Q2 + {e^iXf + ... + Xl^,) ifMet^'). 

The set Z' is the intersection of the sphere S*^ of center and radius - with a 

cylinder based on the extension of Z to R(e). The intersection of Z' with the 
hyperplane Xk+i = is the intersection of Z with the sphere 5*^"^ of center and 

radius -. Denote by tt the projection from R(e)''+^ to R(e)^ 

The following proposition which appears in [22 then relates the connected com- 
ponent of Z with those of Z' and this allows us to reduce the problem of finding 
points on every connected component of a possibly unbounded algebraic set to the 
same problem on bounded algebraic sets. 

Proposition 4.10. Let N he a finite number of points meeting every semi- 
algebraically connected component of Z' . Then n(N) meets every semi-algebraically 
connected component of the extension Ext(Z',R(e)) of Z' to R(e). 



We obtain immediately using Proposition 4.10 a method for finding a point 
in every connected component of an algebraic set. Note that these points have 
coordinates in the extension R(£) rather than in the real closed field R we started 
with. However, the extension from R to R(e) preserves semi-algebraically connected 
components. 

For dealing with possibly singular algebraic sets we define Xi -pseudo-critical 
points of Z(Q, R*-) when Z(Q, R'^) is a bounded algebraic set. These pseudo-critical 
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points arc a finite set of points meeting every semi-algebraically connected compo- 
nent of Z((5,R''"). They are the hmits of the critical points of the projection to the 
Xi coordinate of a bounded nonsingular algebraic hyper-surface defined by a par- 
ticular infinitesimal perturbation, Q, of the polynomial Q. Moreover, the equations 
defining the critical points of the projection on the Xi coordinate on the perturbed 
algebraic set have a very special algebraic structure (they form a Grobner basis [22', 
Section 12.1]), which makes possible efficient computation of these pseudo-critical 
values and points. We refer the reader to [351 Chapter 12] for a full exposition 
including the definition and basic properties of Grobner basis. 

The deformation Q of Q is defined as follows. Suppose that Z((5, R'") is contained 
in the ball of center and radius 1 /c. Let d be an even integer bigger than the 
degree d of Q and let 

(4.2) Gk (d, c) = /(Xf +--- + Xi+Xl + --- + Xl)- {2k ~ 1) , 



(4.3) g = CGfe(J,c) + (l-C)Q. 

The algebraic set Z(Q,R(C) ) is a bounded and non-singular hyper-surface lying 
infinitesimally close to Z(Q,R'^) and the critical points of the projection map onto 
the Xi co-ordinate restricted to Z{Q, R(C) ) form a finite set of points. We take the 



images of these points under lim^ (cf. Definition 4.8 ) and we call the points obtained 
in this manner the Xi-pseudo-critical points of Z((5,R'^). Their projections on the 
Xi-axis are called pseudo-critical values. 

Example 4.11. We illustrate the perturbation mentioned above by a concrete 
example. Let fc = 3 and Q E K[Xi, X2, X3] be defined by 

(4.4) Xl - Xl + X\ + Xl + Xl . 

Then, Z(Q,r3) is a bounded algebraic subset of R'^ shown below in Figure [2] 
Notice that Z(Q,R^) has a singularity at the origin. The surface Z(Q,R^) with 
a small positive real number substituted for C is shown in Figure |3] Notice that 
this surface is non-singular, but has a different homotopy type than Z(Q,R'^) (it 
has three connected components compared to only one of Z((5,R^)). However, the 
semi-algebraic set bounded by Z(Q,R'^) (i.e. the part inside the larger component 
but outside the smaller ones) is homotopy equivalent to Z((5,R^). 



00 



Figure 2. The algebraic set Z((3,R^). 
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Figure 3. The algebraic set Z{Q,R^). 



By computing algebraic representations (see Section 12.4] for the precise 
definition of such a representation) of the pseudo-critical points one obtains for 
any given algebraic set a finite set of points guaranteed to meet every connected 
component of this algebraic set. Using some more arguments from real algebraic 
geometry one can also reduce the problem of computing a finite set of points guar- 
anteed to meet every connected component of the realization of every realizable 
sign condition on a given family of polynomials to finding points on certain alge- 
braic sets defined by the input polynomials (or infinitesimal perturbations of these 
polynomials). The details of this argument can be found in [22, Proposition 13.2]. 

The following theorem which is the best result of this kind appears in jl5j . 

Theorem 4.12. |15j Let Z((5,R'^) be an algebraic set of real dimension k' , where 
Q is a polynomial in Ii[Xi, . . . , Xj,] of degree at most d, and let V C R[^i, . . . , X^] 
be a set of s polynomials with each P CzV also of degree at most d. Let D be the ring 
generated by the coefficients of Q and the polynomials in V. There is an algorithm 
which computes a set of points meeting every semi- algebraically connected compo- 
nent of every realizable sign condition on V over Z{Q, R(e, S)'^). The algorithm has 
complexity 

^k\k - fc') + 1) ^ 4-' ('\ d^C^) = 
j<k' 

in D. There is also an algorithm providing the list of signs of all the polynomials 
of V at each of these points with complexity 



(fc'(fc-fc') + l)s5l4^('')^' 



j<k' 

in D. 



Notice that the combinatorial complexity of the algorithm in Theorem 4.12 de- 
pends on the dimension of the variety rather than that of the ambient space. Since 
we are mostly concentrating on single exponential algorithms in this part of the 
survey, we do not emphasize this aspect too much. 
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4.3. Roadmaps. Theorem |4. 1 2| gives a single exponential time algorithm for test- 
ing if a given semi-algebraic set is empty or not. However, it gives no way of testing 
if any two sample points computed by it belong to the same connected component 
of the given semi- algebraic set, even though the set of sample points is guaranteed 
to meet each such connected component. In order to obtain connectivity infor- 
mation in single exponential time a more sophisticated construction is required - 
namely that of a roadmap of a semi-algebraic set, which is an one dimensional semi- 
algebraic subset of the given semi-algebraic set which is non-empty and connected 
inside each connected component of the given set. Roadmaps were first introduced 
by Canny [30j . but similar constructions were considered as well by Grigoriev and 
Vorobjov [15 and Gournay and Risler [39]. Our exposition below follows that in 
|17l 122] where the most efficient algorithm for computing roadmaps is given. The 
notions of pseudo-critical points and values defined above play a critical role in the 
design of efficient algorithms for computing roadmaps of semi-algebraic sets. 

We first define a roadmap of a semi-algebraic set. We use the following notation. 
We denote by TTi...j the projection, x ^ {xi, . . . , Xj). Given a set 5 C R*^ and 
y &R-' , we denote hy Sy — S (1 T^Z^^jiy)- 

Definition 4.13 (Roadmap of a semi- algebraic set). Let S* C R'^ be a semi- 
algebraic set. A roadmap for 5 is a semi-algebraic set M of dimension at most 
one contained in S which satisfies the following roadmap conditions: 

• RMi For every semi-algebraically connected component D of S, DOM is 
semi-algebraically connected. 

• RM2 For every x £ R and for every semi-algebraically connected component 
D' of S^, D' n M ^ 0. 

We describe the construction of a roadmap RM(Z((5, R'^), A/") for a bounded 
algebraic set Z((5,R'^) which contains a finite set of points Af of Z((5,R'^). A 
precise description of how the construction can be performed algorithmically can 
be found in [35]. We should emphasize here that RM(Z(Q, R*"), A/") denotes the 
semi-algebraic set output by the specific algorithm described below which satisfies 
the properties stated in Definition 4.13 (cf. Proposition 4.14). 

Also, in order to understand the roadmap algorithm it is easier to first concen- 
trate on the case of a bounded and non-singular real algebraic set in R*"' (see Figure 
[4] below). In this case several definitions get simplified. For example, the pseudo- 
critical values defined below are in this case ordinary critical values of the projection 
map on the first co-ordinate. However, one should keep in mind that even if one 
starts with a bounded non-singular algebraic set, the input to the recursive calls 
corresponding to the critical sections (see below) are necessarily singular and thus 
it is not possible to treat the non-singular case independently. 



A key ingredient of the roadmap is the construction of pseudo-critical points and 
values defined above. The construction of the roadmap of an algebraic set con- 
taining a finite number of input points JV of this algebraic set is as follows. We 
first construct A2-pseudo-critical points on Z(Q,R*') in a parametric way along 
the Ai-axis by following continuously, as x varies on the Ai-axis, the A2-pseudo- 
critical points on Z(Q,Ii'^)x. This results in curve segments and their endpoints on 
Z(Q,R'=). The curve segments are continuous semi-algebraic curves parametrized 
by open intervals on the Ai-axis and their endpoints are points of Z((5,R'^) above 




the corresponding endpoints of the open intervals. Since these curves and their 

endpoints include for every .x G R the X2— pseudo-critical points of Z(Q, R'^)^, they 
meet every connected component of Z(<5,R'^)a:- Thus, the set of curve segments 
and their endpoints already satisfy RM2. However, it is clear that this set might 
not be semi-algcbraically connected in a scmi-algcbraically connected component 
and so RMi might not be satisfied. We add additional curve segments to ensure 
connectedness by recursing in certain distinguished hyperplanes defined hy Xi = z 
for distinguished values z. 

The set of distinguished values is the union of the Xi-pseudo-critical values, the 
first coordinates of the input points Af, and the first coordinates of the endpoints 
of the curve segments. A distinguished hyperplane is an hyperplane defined by 
Xi = V, where i; is a distinguished value. The input points, the endpoints of the 
curve segments, and the intersections of the curve segments with the distinguished 
hyperplanes define the set of distinguished points. 

Let the distinguished values he vi < ... < ve- Note that amongst these are the 
Xi-pseudo-critical values. Above each interval {vi,Vi+i) we have constructed a col- 
lection of curve segments Ci meeting every semi-algcbraically connected component 
of Z(g,R'')„ for every v G {vi,Vi-^-l). Above each distinguished value Vi we have a 
set of distinguished points TVj. Each curve segment in Ci has an endpoint in A/i and 
another in Afi+i- Moreover, the union of the J\fi contains J\f. 

We then repeat this construction in each distinguished hyperplane defined by 
Xi = Vi with input Q{vi,X2, . . . , Xk) and the distinguished points in A/^. Thus, we 
construct distinguished values Wi^i, . . . , £(j) of Z(Q(wi, A2, . . . , Afe), R'^~^) (with 
the role of Xi being now played by X2) and the process is iterated until for I = 
(ii, . . . ,ifc_2), 1 < ii < £,...,1 < ik-2 < ({h, ■ ■ ■ ,ik-3), we have distinguished 
values f/j < . . . < vj f(^j^ along the Xk-i axis with corresponding sets of curve 
segments and sets of distinguished points with the required incidences between 
them. 
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Figure 5. The connecting path r(a;) 



The following theorem is proved in [17] (see also [22]). 

Proposition 4.14. The semi- algebraic set RM(Z((5,R'^), A/") obtained by this con- 
struction is a roadmap for Z(Q,I&) containing Af . 

Note that if a; £ Z((3,R''), RM(Z((5, R''), {a:}) contains a path, ^{x), connecting 
a distinguished point p of RM(Z(Q, R*^)) to x. 

4.3.1. The Divergence Property of Connecting Paths. In applications to algorithms 
for computing Betti numbers of semi- algebraic sets it becomes important to examine 
the properties of parametrized paths which are the unions of connecting paths 
starting at a given p and ending at x, where x varies over a certain semi-algebraic 
subset of Z(g,R'=). 

We first note that for any x = (xi, . . . , Xk) S Z((5, R*^) we have by construction 
that RM(Z(g,R'')) is contained in RM(Z(Q, R*"), {a;}). In fact, 

RM(Z(Q, R'-), {x}) = RM(Z(Q, R'^)) U RM(Z(Q, R''%, , J, 

where A4xi consists of (x2, ■ ■ ■ ,Xk) and the finite set of points obtained by inter- 
secting the curves in RM(Z((3, R*"')) parametrized by the Xi-coordinate with the 
hyperplane ttJ" (xi). 

A connecting path 7(x) (with non-self intersecting image) joining a distinguished 
point p of RM(Z(Q,R'')) to x can be extracted from RM(Z(Q, R*-'), {x}). The 
connecting path 7(0;) consists of two consecutive parts, 70(2;) and ri(a;). The path 
7o(a;) is contained in RM(Z((5, R*^)) and the path ri(a;) is contained in Z(Q,R'^)a:j. 
The part 70 (a;) consists of a sequence of sub-paths 7o,Oj ■ ■ • 1 7o,m- Each 70, i is 
a semi-algebraic path parametrized by one of the co-ordinates Xi , . . . , Xk , over 
some interval [ag i, &o,i] with 7o,o(cto,o) = P- The semi-algebraic maps 70,0 j • • • 7 7o,m 
and the end-points of their intervals of definition ag^oi ^o,Oi • • ■ ; c^o,™; ^o,m are all 
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independent of x (up to the discrete choice of the path 7(0;) in RM(Z((5, R*^), {x})), 
except 60, m which depends on xi. 

Moreover, ri(a;) can again be decomposed into two parts 71 (x) and r2(a;) with 
r2(a;) contained in Z(Q, R'^)(^j .j^) and so on. 

If y = (yi, . . . ,yk) G Z((5,R'') is another point such that xi 7^ j/i, then since 
Z{Q,'R'')xi and Z(Q,R'^)j^j^ are disjoint, it is clear that 

RM(Z(g, R'^), {x}) n RM(Z(Q, R^-), {y}) = RM(Z(Q, R'')). 

Now consider a connecting path 7(1/) extracted from RM(Z(Q, R*^), {y}). The im- 
ages of ri(a;) and Ti(y) are disjoint. If the image of 7o(y) (which is contained in 
RM(Z(Q,R'=)) follows the same sequence of curve segments as 70(2;) starting at p 
(i.e. it consists of the same curves segments 70, oi ■ ■ • i7o,m as in 70(2;)), then it is 
clear that the images of the paths 7(x) and 7(y) has the property that they are 
identical up to a point and they are disjoint after it. This is called the divergence 
property in [20,. 

4.3.2. Roadmaps of General Semi-algebraic Sets. Using the same ideas as above and 
some additional techniques for controlling the combinatorial complexity of the al- 
gorithm it is possible to extend the roadmap algorithm to the case of semi-algebraic 
sets. The following theorem appears in |T7l [S^ and gives the most efficient algo- 
rithm for constructing roadmaps. 

Theorem 4.15. Let Q e R[Xi, ...,Xk] with Z(Q,R'') of dimension k' and 

let V C R[^i, . . . , ^fc] &e a set of at most s polynomials for which the degrees of 
the polynomials in V and Q are hounded by d. Let S be a V -semi-algebraic subset 
ofZ{Q,R''). There is an algorithm which computes a roadmap RM(S') for S with 
complexity s'' ) ^j^^^ q generated by the coefficients of Q and the 

elements of V . If D ~ and the bit-sizes of the coefficients of the polynomials 
are bounded by t, then the bit-sizes of the integers appearing in the intermediate 
computations and the output are bounded by rd^''^ ^ . 

Theorem |4.15| immediately implies that there is an algorithm whose output is 
exactly one point in every semi-algebraically connected component of S and whose 
complexity in the ring generated by the coefficients of Q and V is bounded by 
gk +i^o{k )^ jj-^ particular, this algorithm counts the number semi-algebraically 
connected component of S within the same time bound. 



5. Topological Preliminaries 



The purpose of this section is to provide a self-contained introduction to the 
basic mathematical machinery needed later. Some of the topics would be familiar 
to most readers while a few others perhaps less so. The sophisticated reader can 
choose to skip this whole section and proceed directly to the descriptions of the 
various algorithms in the later sections. 

We give a brief review of the concepts from algebraic topology that play a role 
in the results surveyed in this paper. These include the definition of complexes 



of vector spaces (Section 5.2.1 1, definition of cohomology groups of semi-algebraic 



sets (Section 5.3), properties of the Euler-Poincare characteristic of semi-algebraic 
sets 5.3.1), the nerve complex of covers (Section 5.5), a generalization of the nerve 
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complex (Section |5.6[ ), the Mayer- Vietoris double complex and its associated spec- 
tral sequence (Section 5.7), the descent spectral sequence (Section 5.8), and the 
properties of homotopy colimits (Section 5.9). 



5.1. Homology and Cohomology groups. Before we get to the precise defini- 
tions of these groups it is good to have some intuition about them. As noted before 
closed and bounded semi-algebraic sets are finitely triangulable. This means that 
each closed and bounded semi-algebraic set S" C R'^ is homeomorphic (in fact, by 
a semi-algebraic map) to the polyhedron \K\ associated to a finite simplicial com- 
plex K. In fact K can be chosen such that \K\ C R'', and there is an effective 
algorithm (see Theorem 4.5 1 for computing K given S. The simplicial cohomology 
(resp. homology groups) of S are defined in terms of K and are well-defined (i.e 
they are independent of the chosen triangulation which is of course very far from 
being unique). 

Roughly speaking the simplicial homology groups of a finite simplicial complex 
K with coefficients in a field F (which we assume to be Q in this survey) are finite 
dimensional F-vector spaces and measure the connectivity of \K\ in various dimen- 
sions. For example, the zero-th simplicial homology group, IIq(K), has a generator 
corresponding to each connected component of K and its dimension gives the num- 
ber of connected components of \K\. Similarly the first simplicial homology group, 
Hi(i^), is generated by the "one-dimensional holes" of \K\, and its dimension is 
the number of "independent" one-dimensional holes of |iir|. If K is one-dimensional 
(that is a finite graph) the dimension of Hi (A') is the number of independent cycles 
in K. Analogously, the i-th the simplicial homology group, }ii{K), is generated by 
the "z-dimensional holes" of \K\, and its dimension is the number of independent 
i-dimensional holes of \K\. Intuitively an i-dimensional hole is an i-dimensional 
closed surface in K (technically called a cycle) which does not bound any (i + 1)- 
dimensional subset of \K\. 

The simplicial cohomology groups of K are dual (and isomorphic) to the simplicial 
homology groups of K as groups. However, in addition to the group structure 
they also carry a multiplicative structure (the so called cup-product) which makes 
them a finer topological invariant than the homology groups. We are not going 
to use this multiplicative structure. Cohomology groups also have nice but less 
geometric interpretations. Roughly speaking the cohomology groups of K represent 
spaces of globally defined objects satisfying certain local conditions. For example, 
the zero-th cohomology group, H°(if), can be interpreted as the vector space of 
global functions on \K\ which are locally constant. It is easy to see from this 
interpretation that the dimension of H°(if) is the number of connected components 
of K. Similar geometric interpretations can be given for the higher cohomology 
groups, in terms of vector spaces of (globally defined) differential forms satisfying 
certain local condition (of being closed). In literature this cohomology theory is 
referred to as de Rham cohomology theory and it is usually defined for smooth 
manifolds, but it can also be defined for simplicial complexes (see for example [541 
Section 1.3.1]). 



5.1.1. Homology vs Cohomology. It turns out that the cohomological point of view 
gives better intuition in designing algorithms described later in the paper. This is 
our primary reason behind preferring cohomology over homology. Another reason 
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for preferring the cohomology groups over the homology groups is that their inter- 
pretations continue to make sense in applications outside of semi-algebraic geometry 
where the notions of holes is meaningless (for instance, think of algebraic varieties 
defined over fields of positive characteristics) but the notion of global functions (or 
for instance differential forms) continue to make sense. 

5.2. Definition of the Cohomology Groups of a Simplicial Complex. We 

now give precise definitions of the cohomology groups of simplicial complexes. 

In order to do so we first need to introduce some amount of algebraic machinery 
- namely the concept of complexes of vector spaces and homomorphisms between 
them. 

5.2.1. Complex of Vector Spaces. A complex of vector spaces is just a sequence of 
vector spaces and linear transformations satisfying the property that the composi- 
tion of two successive linear transformations is 0. 
More precisely 

Definition 5.1 (Complex of Vector Spaces). A sequence {C}, p G Z, of Q- vector 
spaces together with a sequence {S^} of homomorphisms (5^* : — *■ G^^^ (called 
differentials) for which 

(5.1) 6P+^ oSP^O 
for all p is called a complex. 

The most important example for us of a complex of vector spaces is the co-chain 
complex of a simplicial complex K denoted by C*{K). It is defined as follows. 

Definition 5.2 (Simplicial cochain complex). For each p > 0, GP{K) is a linear 
functional on the Q-vector-space generated by the p-simplices of K. Given (j) € 
CP{K), 6P{(f)) is specified by its values on the {p + l)-dimensional simplices of K. 
Given a, (p + l)-dimensional simplex a — [oq, . . . , flp+i] of K 

P+i 

(5.2) {SP(l)){[ao, ap+i]) = 'V(-l) V([ao, . . . , d,, . . . , Cp+i]), 



i=0 



where " denotes omission. 



5.1 



Notice that each [gq, . . . ,di, . . . , a^+i] is a p-dimensional simplex of K and since 
(j) G CP{K), 0([ao, . . . ,di, . . . , Op+i]) S Q is well-defined. It is an exercise now to 
check that the homomorphisms 6^ : Gp{K) Gp^^{K) indeed satisfy Eqn. 
the definition of a complex. 

Now let be a simplicial complex and L C K a. sub-complex of K - we will 
denote such a pair simply by {K, L). Then for each p > we have that Cp(L) C 
GP{K) and we denote by Cp{K,L) the quotient space Cp{K)/Cp{L). It is now an 
easy exercise to verify that the differentials Sp in the complex Cp{K) descend to 
CP{K, L) and we define 

Definition 5.3 (Simplicial cochain complex of a pair). The simplicial cochain 
complex of the pair {K, L) to be the complex C'{K, L) whose terms, Gp{K, L), and 
differentials, Sp, are defined as above. 



Often, particularly in the context of algorithmic applications it is more econom- 
ical to use cellular complexes instead of simplicial complexes. We recall here the 
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definition of a finite regular cell complex referring the reader to standard sources 
in algebraic topology for more in-depth study of cellular theory (see |69l pp. 81]). 

Definition 5.4 (Regular cell complex). An ^^-dimensional cell in R*^ is a subset of 
R'° homeomorphic to 5^(0, 1). A regular cell complex S in R*^ is a finite collection 
of cells satisfying the following properties: 

(1) If ci, C2 e S, then either ci n C2 = or ci C dc2 or C2 C dci. 

(2) The boundary of each cell of E is a union of cells of E. 

We denote by |E| the set [J c. 

Remark 5.5. Notice that every simplicial complex K may be considered as a regular 
cell complex whose cells are the closures of the simplices of K. 

As in the case of simplicial complexes it is possible to associate a complex, C*(E) 
(the co-chain complex of K) , to each regular cell complex K which is defined in an 
analogous manner. In order to avoid technicalities we omit the precise definition 
of this complex referring the interested reader to [69l pp. 82] instead. We remark 
that the dimension of 0^(2) is equal to the number of p-dimensional cells in S and 
the matrix entries for the differentials in the complex with respect to the standard 
basis comes from {0, 1, —1} just as in the case of simplicial co-chain complexes. 

The advantage of using cell complexes instead of simplicial complexes can be 
seen in the following example. 

Example 5.6. Consider the unit sphere S*^ C R'^^^ For < j < fc -t- 1 and 

£ e {+, -} let 

(5.3) = {x e S'^ I Xo = --- = Xj_i = 0, eXj > 0}. 

Then it is easy to check that each is a fc — j dimensional cell and the collection, 

— I < j < fc,£ e {^,+}} is a regular cell complex with \Y.k\ — S*"' (see 
Figure Ig] for the case k = 2). 




Figure 6. Cell decomposition of 
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Notice that #5]^: = 2k. On the other hand if we consider the sphere as home- 
omorphic to the boundary of a standard (k + l)-dimensional simplex, then the 
corresponding simplicial complex will contain (2''+^ — 2) simplices (which is expo- 
nentially large in fc). 

We now associate to each complex, C*, a sequence of vector spaces, W{G'), 
called the cohomology groups of C*. Note that it follows from Eqn. |5.1| that for a 
complex C* with differentials Sp : CP ^ Cp+^ the subspace BP{C') = lin{6P-^) C 
CP is contained in the subspace ZP{C') = Ker{6P) C C^. The subspaces BP{C') 
(resp. ZP{C*)) are usually referred to as the co-boundaries (resp. co-cycles) of the 
complex C*. Moreover, 

Definition 5.7 (Cohomology groups of a complex). The cohomology groups, H''(C*) , 
are defined by 

(5.4) }iP{C') = ZP{C')/BP{C'). 

We will denote by H*(C*) the graded vector space ®pHP(C*). 

Note that the cohomology groups, W{C*), are all Q- vector spaces (finite dimen- 
sional if the vector spaces C^'s are themselves finite dimensional). 

Definition 5.8 (Exact sequence). A complex C* is called acyclic and the corre- 
sponding sequence of vector space homomorphisms is called an exact sequence if 
H*(C*) = 0. 

Applying Definition |5.7| to the particular case of the co-chain complex of a sim- 



plicial complex K (cf. Definition 5.2) we obtain 



5.2.2. Cohomology of a Simplicial Complex. 

Definition 5.9 (Cohomology of a simplicial complex). The cohomology groups of 
a simplicial complex K are by definition the cohomology groups, 'HP{C'{K)), of its 
co-chain complex. 

Similarly, given a pair of simplicial complexes {K, L) , we define 

Definition 5.10 (Cohomology of a pair). The cohomology groups of the pair {K, L) 
are by definition the cohomology groups, W{C*{K,L)), of its co-chain complex. 

Example 5.11. Let A„ be the simplicial complex corresponding to an n-simplex. 
In other words the simplices of A„ consist of [io, . . . ,ie\,0 < iq < ■ ■ ■ < ie < n. The 



polyhedron |A„| is just the n-dimensional simplex. Then using Definition 5.9 
can verify that 

ff (A„) = Q, i = 0, 
ff (A„) = 0, I > 0. 

Example 5.12. Let 9A„ be the simplicial complex corresponding to the boundary 
of the n-simplex. In other words the simplices of 9A„ consist of [ig, . . . , i^], < iq < 
■ ■ ■ < ig < n,£ < n. Then again by a direct application of Definition |5.9| one can 
verify that 

H'(aA„) = Q, i = 0,n- 1 
W{dAn) = 0, else. 
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The above examples serve to confirm our geometric intuition behind the homol- 
ogy groups of the spaces |A„| and |i9A„| explained in Section 5.1 above - namely 
that they are both connected and |A„| has no holes in dimension > 0, and |i9A„| 
has a single (n — l)-dimensional hole. 

Example 5.13. It is also an useful exercise to verify that 

ff(A„,aA„) =Q, i = Q,n 

ff(A„,aA„)=0, else. 



Remark 5.14. Example |5.13| illustrates that for "nice spaces" of the kind we consider 
in this paper (such as regular cell complexes) the cohomology groups of a pair 
{K, L) are isomorphic to the cohomology groups of the quotient space |if|/|L|. For 
instance, the above example illustrates the fact that the topological quotient of an 
n-dimensional ball by its boundary is the n-dimensional sphere. 

5.2.3. Homomorphisms of Complexes. We will also need the notion of homomor- 
phisms of complexes which generalizes the notion of ordinary vector space homo- 
morphisms. 

Definition 5.15 (Homomorphisms of complexes). Given two complexes, C* = 
(C^,(SP) and D* = (D^,J^), a homomorphism of complexes, (f>' : C ^ D*, is a 
sequence of homomorphisms : ^ T)^ for which (5^ o 0^ = (pP'^^ o 5^ for all p. 
In other words the following diagram is commutative. 



5" 



4. Qp+l 



(5.5) 



A homomorphism of complexes (/>* : C* — > D* induces homomorphisms (j/ : 
H*(C*) — > H'(D*) and we will denote the corresponding homomorphism between 
the graded vector spaces H*(C*), H*(D*) by 0*. 



Definition 5.16 (Quasi-isomorphism). The homomorphism 
isomorphism if the homomorphism <j)* is an isomorphism. 



is called a quasi- 



Having introduced the algebraic machinery of complexes of vector spaces, we now 
define the cohomology groups of semi-algebraic sets in terms of their triangulations 
and their associated simplicial complexes. 

5.3. Cohomology Groups of Semi-algebraic Sets. A closed and bounded semi- 
algebraic set 5' C R*^ is semi-algebraically triangulable (see Theorem 



4.5 



above) . 



Definition 5.17 (Cohomology groups of closed and bounded semi-algebraic sets). 
Given a triangulation, h : \K\ — > S, where K is a, simplicial complex, we define 
the i-th simplicial cohomology group of S, by H*(5') = W{C*{K)), where C*{K) is 
the co-chain complex of K. The groups H*(S') are invariant under semi-algebraic 
homeomorphisms (and they coincide with the corresponding singular cohomology 
groups when R = M). We denote by bi{S) the «-th Betti number of S (i.e. the 
dimension of H'(S') as a vector space). 



28 



SAUGATA BASU 



Remark 5.18. For a closed but not necessarily bounded semi-algebraic set C R*^ 
we will denote by H*(S') the i-th simplicial cohoniology group of S* n Bk{0,r) for 
sufficiently large r > 0. The sets S D Bk{0,r) are semi-algebraically homeomorphic 
for all sufficiently large r > and hence this definition makes sense. (The last 
property is usually referred to as the local conic structure at infinity of semi-algebraic 
sets [22l Theorem 5.48]). The definition of cohomology groups of arbitrary semi- 
algebraic sets in R requires some care and several possibilities exist and we refer 
the reader to Section 6.3] where one such definition is given which agrees with 
singular cohomology in case R = M. 

5.3.1. The Euler-Poincare Characteristic: Definition and Basic Properties. An use- 
ful topological invariant of semi-algebraic sets which is often easier to compute than 
their Betti numbers is the Euler-Poincare characteristic. 

Definition 5.19 (Euler-Poincare characteristic of a closed and bounded semi-al- 
gebraic set). Let S C R*", be a closed and bounded semi-algebraic set. Then the 
Euler-Poincare characteristic of S is defined by 

(5.6) x{s)=j2(-^yb^(sy 

i>0 

From the point of view of designing algorithms, it is useful to define Euler- 
Poincare characteristic also for locally closed semi-algebraic sets. A semi-algebraic 
set is locally closed if it is the intersection of a closed semi-algebraic set with an open 
one. A standard example of a locally closed semi-algebraic set is the realization, 
TZ{cr), of a sign-condition cr on a family of polynomials. 

We now define Euler-Poincare characteristic for locally closed semi-algebraic sets 
in terms of the Borel-Moore cohomology groups of such sets (defined below). This 
definition agrees with the definition of Euler-Poincare characteristic stated above 
for closed and bounded semi-algebraic sets. They may be distinct for semi-algebraic 
sets which are closed but not bounded. 

Definition 5.20. The simplicial cohomology groups of a pair of closed and bounded 
semi-algebraic sets T C S" C R*^ are defined as follows. Such a pair of closed and 



bounded semi-algebraic sets can be triangulated (cf. Theorem 4.5 1 using a pair of 
simplicial complexes {K, A) where A is a sub-complex of K. The p-th simplicial 
cohomology group of the pair {S,T), }iP{S,T), is by definition to be IiP{K,A). 
The dimension of W{S,T) as a Q- vector space is called the p-th Betti number of 
the pair {S, T) and denoted 6p(S', T). The Euler-Poincare characteristic of the pair 
{S,T) is 

X(5,T) = ^(-1)^5.(5, T). 

i 

Definition 5.21 (Borel-Moore cohomology group). The p-th Borel-Moore coho- 
mology group of 5* C R*^, denoted H^j^^(5'), is defined in terms of the cohomology 
groups of a pair of closed and bounded semi-algebraic sets as follows. For any r > 
let Sr = S'ni?fe(0, r). Note that for a locally closed semi-algebraic set S both Sr and 
Sr\Sr are closed and bounded, and hence W{SrTSr\Sr) is well defined. Moreover, 
it is a consequence of the local conic structure at infinity of semi-algebraic sets (see 
above) that the cohomology group HP{Sr, Sr \ Sr) is invariant for all 



5.18 



Remark 

sufficiently large r > 0. We define H^^,^(S') = HP(S^,SV \ Sr) for r > sufficiently 
large and it follows from the above remark that it is well defined. 
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The Borel-Moore cohomology groups are invariant under semi-algebraic homeo- 
morphisms (see [2^ . It also follows clearly from the definition that for a closed and 
bounded semi-algebraic set the Borel-Moore cohomology groups coincide with the 
simplicial cohomology groups. 

Definition 5.22 (Borel-Moore Euler-Poincare characteristic). For a locally closed 
semi-algebraic set S we define the Borel-Moore Euler-Poincare characteristic by 

(5.7) ^BM^s-^^j2{-irbr'{s) 

1=0 

where bf^^{S) denotes the dimension of H^^,j(S'). 

Since the Borel-Moore Euler-Poincare characteristic might not be very familiar, 
the reader is encouraged to compute it in a few simple examples. In particular, one 
should check that for the half-open interval [0, 1) which is locally closed we have 

(5.8) x'''''([0,l)) = 0. 
We also have 

(5.9) x''*'((0,l))--l, 
and more generally, 

(5.10) x^'^(Bfc(0,l)) = (-l)^ 
If S is closed and bounded then x^^^{S) = xi^)- 

The Borel-Moore Euler-Poincare characteristic has the following additivity prop- 
erty (reminiscent of the similar property of volumes) which makes them particularly 
useful in algorithmic applications (see for example Section 7.3 below). 

Proposition 5.23. Let X, Xi and X2 be locally closed semi- algebraic sets such 
that 

Xi u X2 = X, Xi n = 0. 

Then 

(5.11) x^""(^) = X^""(^i) + X^'-'(^2). 

Since for closed and bounded semi-algebraic sets, the Borel-Moore Euler-Poincare 
characteristic agrees with the ordinary Euler-Poincare characteristic, it is easy to de- 
rive the following additivity property of the Euler-Poincare characteristic of closed 
and bounded sets. 

Proposition 5.24. Let Xi and X2 be closed and bounded semi- algebraic sets. Then 

(5.12) x{Xi U X2) = x{Xi) + x{X2) - x{Xi n X2). 

Note that Proposition |5.24| is an immediate consequence of Proposition |5 . 23 1 once 
we notice that the sets Yi — Xi\ {Xi n X2) and I2 = -^^2 \ (-'^i H X2) are locally 
closed, the set Xi U X2 is the disjoint union of the locally closed sets Yi,F2 and 
XinX2, and 

X^'\Y,) = x^^(X,) - x^^(Xi n X2) = x{X,) - x{Xi n X^), for i = l,2. 

More generally by applying Proposition |5.24| inductively we get the following 
inclusion-exclusion property of the (ordinary) Euler-Poincare characteristic. 
For any n £ Z>o we denote by [n] the set {1, . . . , n}. 
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Proposition 5.25. Let Xi,...,Xn be closed and bounded semi-algebraic sets. 
Then denoting by Xj the semi- algebraic set Xi for / C [n], we have 

iei 

(5.13) x( U = E (-1)'*'+'^ xiXi). 

ie[n] IC[n] 

5.4. Homotopy Invariance. The cohomology groups of semi-algebraic sets as 



defined above (Definition 5.171 are obviously invariant under semi-algebraic home- 
omorphisms. But, in fact, they are invariant under a weaker equivalence relation 
- namely, semi-algebraic homotopy equivalence (defined below). This property is 
crucial in the design of efficient algorithms for computing Betti numbers of semi- 
algebraic sets since it allows us to replace a given set by one that is better behaved 
from the algorithmic point of view but having the same homotopy type as the orig- 
inal set. This technique is ubiquitous in algorithmic semi-algebraic geometry and 
we will see some version of it in almost every algorithm described in the following 



sections (cf. Example 4.11) 



Remark 5.26. The reason behind insisting on the prefix "semi-algebraic" with re- 
gard to homeomorphisms and homotopy equivalences here and in the rest of the 
paper, is that for general real closed fields, the ordinary Euclidean topology could 
be rather strange. For example, the real closed field, Maig, of real algebraic numbers 
is totally disconnected as a topological space under the Euclidean topology. On the 
other hand, if the ground field R = M, then we can safely drop the prefix "semi- 
algebraic" in the statements made above. However, even if we start with R = M, in 
many applications described below we enlarge the field by taking non-archimedean 



extensions of R (see Section 4.2.1 ), and the remarks made above would again apply 
to these field extensions. 

Definition 5.27 (Semi-algebraic homotopy). Let X, Y be two closed and bounded 
semi-algebraic sets. Two semi-algebraic continuous functions f,g : X Y are 
semi- algebraically homotopic, f '^sa 9, if there is a continuous semi-algebraic func- 
tion F : X X [0, 1] ^ r such that F{x, 0) = f{x) and F{x, 1) = g{x) for all x€ X. 

Clearly, semi-algebraic homotopy is an equivalence relation among semi-algebraic 
continuous maps from X to Y . 

Definition 5.28 (Semi-algebraic homotopy equivalence). The sets X,Y are semi- 
algebraically homotopy equivalent if there exist semi-algebraic continuous functions 
f : X ^Y, g -.Y ^ X such that go f r^^^ Wx, f o 9 ^sa Wy. 

We have 

Proposition 5.29 (Homotopy Invariance of the Cohomology Groups). Let X,Y 

be two closed and bounded semi-algebraic sets of R*^ that are semi- algebraically 
homotopy equivalent. Then, H*(X) = Il*(Y). 

5.5. The Leray Property and the Nerve Lemma. It clear from the definition 



of the cohomology groups of closed and bounded semi-algebraic sets (Definition 5.17 
above) the Betti numbers of such a set can be computed using elementary linear 
algebra once we have a triangulation of the set. However, as we have seen before 



(cf. Theorem 4.51, triangulations of semi-algebraic sets are expensive to compute, 



requiring double exponential time. 
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One basic idea that underlies some of the recent progress in designing algo- 
rithms for computing the Betti numbers of semi-algebraic sets is that the coho- 
mology groups of a semi-algebraic set can often be computed from a sufficiently 
well-behaved covering of the set without having to triangulate the set. 

The idea of computing cohomology from "good" covers is an old one in algebraic 
topology and the first result in this direction is often called the "Nerve Lemma" . In 
this section we give a brief introduction to the Nerve Lemma and its generalizations. 

We first define formally the notion of a cover of a closed, bounded semi-algebraic 
set. 

Definition 5.30 (Cover). Let S C R*" be a closed and bounded semi-algebraic set. 
A cover, C{S), of S consists of an ordered index set, which by a slight abuse of 
language we also denote by C{S), and a map that associates to each a £ C{S) a 
closed and bounded semi-algebraic subset Sa C S such that 

S = Sa- 

aeC{S) 

Remark 5.31. Even though the notation for a cover might seem unnecessarily heavy 
at the moment it will prove useful later on the paper when we discuss non-Leray 



covers (see Section 5.6 below) 



For ckg, . . . , ap, e C(S'), we associate to the formal product, ag • • • Op, the closed 
and bounded semi-algebraic set 



(5.14) Sao-a,^ Saor\---nSa,. 

Recall that the 0-th simplicial cohomology group of a closed and bounded semi- 
algebraic set X, H°(X), can be identified with the Q- vector space of Q- valued 
locally constant functions on X. Clearly the dimension of H°(X) is equal to the 
number of connected components of X. 

For ao, ai, . . . , ap, (3 £ C{S), and (3 ^ {ao, . . . , otp}, let 

^Qo,...,ap;/3 • H°(S'Qf,...Q,p) > H°(S'QQ...ap.^) 

be the homomorphism defined as follows. Given a locally constant function, e 
,;/3(0) is the locally constant function on S'ao - ap /3 obtained by 
restricting 4) to Sao--aj,-i3- 

We define the generalized restriction homomorphisms 

■ H'^(S'qq...Qp) > H°(S'c,o...Qp^j) 

by 

(5.15) <^^(</')ao---c<p+l ~ ^ ^ (~1) ^Qo---c<i---c«p+i;ai (0ao---cfi---ap4,i )j 

0<i<p+l 

where (f) E ®ao< - <ap<£C{S) 

H°(^o„...„J and fao---di---aj,+i;ai is the restriction ho- 
momorphism defined previously. The sequence of homomorphisms 6^ gives rise to 
a complex, L*{C{S)), defined by 



(5.16) 



LP{C{S)) = h0(^„„...„J, 

ao<---<ap,ai£C{S) 
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with the differentials : LP{C{S)) LP+'^{C{S)) defined as in Eqn. (|5.15[) 



Definition 5.32 (Nerve complex). The complex L*{C{S)) is called the nerve com- 
plex of the cover C{S). 

For £ > we will denote by L*{C{S)) the truncated complex defined by 

mC{S))^LP{C{S)), 0<p<i, 
^0, p>i. 

Notice that once we have a cover of S and we identify the connected components 
of the various intersections, Sa„---ap, we have natural bases for the vector spaces 

LP{CiS)) = hO(5„„...„J 

ao<---<ap,ai£C{S) 

appearing as terms of the nerve complex. Moreover, the matrices corresponding 
to the homomorphisms 6^ in this basis depend only on the inclusion relationships 
between the connected components of 

Sao---ap^i and those of ScxQ---ap- 

Definition 5.33 (Leray Property). We say that the cover C{S) satisfies the Leray 
property if each non-empty intersection Sao---ap is contractible. 

Clearly, in this case 

^ 0, if = 0. 

It is a classical fact (usually referred to as the Nerve Lemma) that 

Theorem 5.34 (Nerve Lemma). Suppose that the cover C{S) satisfies the Leray 
property. Then for each i > 0, 

ff(L*(C(5))) ^^ff(5). 

(See for instance [BT for a proof.) 



Remark 5.35. There are several interesting extensions of Theorem 5.34 (Nerve 
Lemma). For instance, if the Leray property is weakened to say that each t-ary 
intersection is {k — t+ l)-connected, then one can conclude that the nerve complex 
is fc-connected. We refer the reader to the article by Bjorner [55] for more details. 

Notice that Theorem |5 . 34| gives a method for computing the Betti numbers of S 
using linear algebra from a cover of S by contractible sets for which all non-empty 
intersections are also contractible, once we are able to test emptiness of the various 
intersections Sao---ap- 

Now suppose that each individual member, Sa„ , of the cover is contractible, but 
the various intersections Sao---ap are not necessarily contractible for p > I. Theorem 



5.34 does not hold in this case. However, the following theorem is proved in |20] 
and underlies the single exponential algorithm for computing the first Betti number 
of semi-algebraic sets described there. 

Theorem 5.36. [20] Suppose that each individual member, Sag, of the cover C{S) 
is contractible. Then, 

ff(L*(C(5)))-ff(5), fori = Q,l. 
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Remark 5.37. Notice that from a cover by contractible sets Theorem |5.36| allows 
us to compute using linear algebra, bQ^S) and bi(S), once we have identified the 
non-empty connected components of the pair-wise and triple-wise intersections of 
the sets in the cover and their inclusion relationships. 



Example 5.38. We illustrate Remark 5.37 with a simple example. 

Consider the following set S depicted in Figure [7]below and let C{S) = {0, 1,2} 
and the corresponding sets 5*0, ^i, 5*2 are the three edges as shown in in the figure. 




Figure 7. Example illustrating Theorem 5.36 



Notice that each pair-wise and triple-wise intersections in this case has two con- 
nected components. Let us construct the complex L2(C(5)). We have 

(5.17) L°(C(^))= HO(5o)®HO(5i)®HO(5i) ?^^Q®Q®Q, 

aoeC(S) 

(5.18) Li(C(5)) = H°(S'oi)eH°(5o2)©H"(^i2) = QeQeQeQeQeQ, 

ao<aieC{S) 

and 

(5.19) l2(C(^))= HO(S'oi2) -Q®Q. 

ao<ai <a2SC(S) 

We now display the matrices Mq and Mi corresponding to the homomorphisms 
(5° and respectively (with respect to the obvious basis corresponding to the 
connected components of the various intersections). 

We have 



(5.20) 

and 
(5.21) 



Mo 



Ml 





/ 1 


-1 





\ 




1 


-1 









1 





-1 






1 





-1 









1 


-1 








1 


-1 


/ 


1 





-1 





1 





1 





-1 
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It is now easy to verify that rank(Mo) — 2 and rank(Mi) = 2. We derive 
applying Theorem |5.36| that 

bo{S) = dimkerMo 

= dim(L°(C(S'))) - rank(Mo) 
= 3-2 
= 1, 

and 

bi{S) = dimker(Mi) - dimlm(Mo) 

= dini(L^ (C(S'))) - rank(Mi) - rank(Mo) 

=6-2-2 

= 2. 

We refer the reader to [33] for more comphcated higher dimensional examples of a 
similar nature. 

Remark 5.39. It is easy to see that if we extend the complex in Theorem |5.36| 
by one more term, that is consider the complex, 1j*{C{S)), then the cohomology 
of the complex does not yield information about 11^(5). Just consider the cover 
of the standard sphere C R"^ and the cover {Hi,H2} of where Hi,H2 are 
closed hemispheres meeting at the equator. The corresponding complex, Ij*(C), is 
as follows. 

^ H"(Hi) 0HO(i/2) ^ H"(i/i n i/2) ^ ^ 

Clearly, H^(L3(C)) 9^ H^(S^), and indeed it is impossible to compute bi{S) just 
from the information on the number of connected components of intersections of 
the sets of a cover of S by contractible sets for i > 2. For example the nerve 
complex corresponding to the cover of the sphere by two hemispheres is isomorphic 
to the nerve complex of a cover of the unit segment [0, 1] by the subsets [0,1/2] and 
[1/2, 1], but clearly lP{S^) = Q, while I1^{[0, 1]) = 0. 

5.6. Non-Leray Covers. In the design of algorithms for computing covers of semi- 
algebraic sets it is often difficult to satisfy the full Leray property. In order to utilize 
covers not satisfying the Leray property it is necessary to consider a generalization 
of the nerve complex. However, before we can describe this generalization we need 
to expand slightly the algebraic machinery at our disposal. 

We first introduce the notion of a double complex which is in essence a complex 
of complexes. 

Definition 5.40 (Double complex). A double complex is a bi-graded vector space 

with co-boundary operators d : CP-? CP'«+i and 5 : CP'^ ^ Cp+1'« and such 
that doS + Sod — (see diagram below). We say that C*'* is a first quadrant 
double complex if it additionally satisfies the condition that C^ '^ = if either p < 
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or q < 0. Double complexes lying in other quadrants are defined in an analogous 
manner. 



q0,2 q1,2 q2,2 



qO,1 Ql.l q2,1 



Definition 5.41 (The Associated Total Complex). The complex defined by 

Tot"(C*'*)= Cf'«, 

p+q=7l 

with differential 

D" = d+i-lfd : Tot"(C*'*) — > Tot"+\C*'*), 

p+q=n 

is denoted by Tot*(C*'*) and called the associated total complex of C*'*. 

Associated to double complexes, or more accurately to their filtrations, is another 
algebraic object that is quite ubiquitous in modern algebraic topology. 

Definition 5.42 (Spectral Sequence). A (cohomology) spectral sequence is a se- 
quence of bi-graded (this is a direct sum of vector subspaces indexed by Z x Z) com- 
plexes {El'^ \ i,j,r € Z,r > a} endowed with differentials dp^ : Ep^ j^i+r,j-r+i 
such that (dr)^ = for all r. Moreover, we require the existence of isomorphism 
between the complex -Br+i and the homology of E,. with respect to d^: 

ker d^^-^ 

El'l, = HdAK^n = ^ 

'■+1 ^^'■'^-'■+1 (ei+'-'^-'-+^) 

The spectral sequence is called a first quadrant spectral sequence (see Figure [S]) 
if the initial complex Ea lies in the first quadrant, i.e. i?*'^ = whenever ij < 0. 
In that case, all subsequent complexes Er also lie in the first quadrant. Since the 
differential d*'-' maps outside of the first quadrant for r > i, the homomorphisms of 
a first quadrant spectral sequence dr are eventually zero, and thus the groups Ep^ 
are all isomorphic to a fixed group E^^ for r large enough, and we say the spectral 
sequence is convergent. 
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1 k 




p+q=e p+q=e+l 



Figure 8. dr : EP'I i;p+r,9-r-+i 

There are two spectral sequences, 'E^'"^ ,"E^''' , (corresponding to taking row- wise 
or column-wise filtrations respectively) associated with a first quadrant double com- 
plex C*'* which will be important for us. Both of these converge to H*(Tot*(C*'*)). 
This means that the homomorphisms, dr are eventually zero, and hence the spectral 
sequences stabilize and 

(5.22) - "E^J - (Tof(C-'-)) 

p+q=i p+q=i 

for each z > 0. 

The first terms of these are 

(5.23) 'El = Rd{C'''),'E2 = H5Hd(C*'*), 

and 

(5.24) "El = Rs{C'n,"E2 = ^d^siC'')- 

Given two (first quadrant) double complexes, C*'* and C*'*, a homomorphism 
of double complexes, 

. (J.,. — ^ 

is a collection of homomorphisms (j)^''^ : G^''^ — > C^'*^ such that all the following 
diagrams commute. 
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QP,q ; (^p,q+l 

A homomorphism of double complexes 

(/)'•' : C*'* C*'*, 

induces an homomorphism of the corresponding total complexes which we will 
denote by 

Tof(0*^*) : Tof(C*'*) — > Tof(C*^*). 

It also induces homomorphisms '(ps ■ 'Eg — > 'Eg (respectively, "4>s : "Eg — > "Eg) 
between the associated spectral sequences (corresponding either to the row-wise or 
column- wise filtrations). For the precise definition of homomorphisms of spectral 
sequences see [52 . We will need the following useful fact (see [52, pp. 66]). 

Theorem 5.43 (Comparison Theorem). If'(j)g (respectively, "4>g) is an isomor- 
phism for some s > 1 then 'E^.''' and 'E^''^ (respectively, " E^''^ and "E'^'"^) are 
isomorphic for all r > s. In particular the induced homomorphism 

Tof(<?!>*'*) : Tof(C*^*) — > Tof(C*'*) 

is a quasi-isomorphism. 

Having introduced the definitions of double complexes and spectral sequences 
above, we now describe two particular double complexes that are of interest to us 
in this paper. 

5.7. The Mayer- Vietor is Double Complex and its Associated Spectral 
Sequence. Let Ai, . . . , An be sub-complexes of a finite simplicial complex A such 
that A — Ai U ■ ■ ■ U An- Note that the intersections of any number of the sub- 
complexes, Ai, is again a sub-complex of A. We denote by A^g.-.a^ the sub-complex 
Aao n • ■ ■ n Aa^j ■ 

Definition 5.44 (The Generalized Mayer- Vietoris Exact Sequence). The general- 
ized Mayer- Vietoris sequence is the following exact sequence of vector spaces. 

> C*(^) > G'{Aag) > G'{Aao.ai) 

^3 C'{Aao...ap) > ^ C {Aag...ap+i) 

l<C(o<---<Qp<n l<ao<---<ap+i<n 

where r* is induced by restriction and the homomorphisms 5^'* are defined as 
follows. 

Given an w S 0qo<...<qj, C'^{^ao---ap) we define SP''^{u!) as follows: 
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First note that 6^''^uj e C'(^c((,...Qp^j), and it suffices to define 

ao<---<ap+l 

for each (p + 2)-tuplc f < ao < • • • < ctp+i ^ Note that, {5^'''Lo)ao,...,ap+i is 
a hnear form on the vector space, Cq{Aci^^...aj,j^-^), and hence is determined by its 
values on the g-simphces in the complex ^^^...0,^^^ . Furthermore, each g-simplex, 
s e ^ao -ap+i is automatically a simplex of the complexes 



A 



1, < i <p+ 1. 



We define 



0<j<p+l 

The fact that the generalized Mayer- Vietoris sequence is exact is classical (see 
[5T] or [7] for example). 

We now define the Mayer- Vietoris double complex of the complex A with respect 
to the sub-complexes Aa^,! < ao < n, which we will denote by N'''{A) (we 
suppress the dependence of the complex on sub-complexes Aa^ in the notation 
since this dependence will be clear from context). 

Definition 5.45 (Mayer- Vietoris Double Complex). The Mayer- Vietoris double 
complex of a simplicial complex A with respect to the sub-complexes Aa^ , 1 < 
ao < n. A/"*'* (A), is the double complex defined by 

l<Qo<-<ap<n 

The horizontal differentials are as defined above. The vertical differentials are those 
induced by the ones in the different complexes, G'{Aaa---a^)- 
N'''{A) is depicted in the following figure. 



(5.25) 



ao<Q:i 



ao<Q:i 



ao<Q:i 
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Remark 5.46. There is also dual version of the Mayer- Vietoris double complex where 
the unions and intersections are inter-changed and the directions of the arrows get 
reversed. The reader is referred to [T for more detail. 

Finally, for complexity reasons it is often useful to consider truncations of the 
Mayer- Vietoris double complex. 

Definition 5.47 (Truncated Mayer- Vietoris Double Complex). For any i > 0, we 
denote by M'''{A) the following truncated complex. 

A/7'«(^) = AAP'9(A), Q<v + q<t, 
Nf''^{A) = 0, otherwise. 

The following proposition is classical (see [61] or [7] for a proof) and follows from 
the exactness of the generalized Mayer- Vietoris sequence. 

Proposition 5.48. The spectral sequences, 'E^/'Ej., associated to J\f*'*{A) con- 
verge to H*(^) and thus, 

H*(Tof(AA*'*(yl))) = H*(A). 

Moreover, the homomorphism 

tP' : C'{A) ^ Tot'{Af'''{A)) 

induced by the homomorphism r* (in the generalized Mayer- Vietoris sequence) is a 
quasi-isomorphism. 

We denote by C*^^{A) the truncation of the complex C*{A) after the {£-\- l)-st 
term. As an immediate corollary we have that, 

Corollary 5.49. For any £ > 0, the homomorphism 

(5.26) r^+i ■■ C,Vi(A) ^ Tot'iJV;^\iA)) 

induced by the homomorphism r* (in the generalized Mayer- Vietoris sequence) 
is a quasi-isomorphism. Hence, for < i < i, 

ff(Tof(AA;_,-(^)))=ff(^). 

Remark 5.50. Notice that in the truncated Mayer- Vietoris double complex, Af*''{A), 
the 0-th column is a complex having at most t-\-l non-zero terms, the first column 
can have at most t non-zero terms, and in general the z-th column has at most 
t-\-l—i non-zero terms. This observation along with the fact that, each term in the 
double complex Af'''{A) depends on tuples of at most t -\- 1 of the A^ 's at a time, 
play a crucial role in the inductive arguments used in the design of single expo- 
nential time algorithm for computing the first few Betti numbers of semi-algebraic 
sets. 

5.8. The Descent Double Complex and its Associated Spectral Sequence. 

For the algorithmic problem of computing the Betti numbers of projections of semi- 
algebraic sets, another spectral sequence plays an important role. 

Definition 5.51 (Locally split maps). A continuous surjection / : A" — > F is called 
locally split if there exists an open covering U of Y such that for all U G U, there 
exists a continuous section a : U X of f , i.e. cr is a continuous map such that 
/(o-(y)) = y for all yeU. 
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In particular, if X is an open semi-algebraic set and f : X Y is a, projection, 
the map / is obviously locally split. For any semi-algebraic surjection f : X ^ Y, 
we denote by W^{X) the {p + l)-fold fibered power of X over /, 



W^{X) = {{xo, ...,xp)e XP+' I f{xo) = • • • = f{xp)}. 



The map / induces for each p > 0, a map from Wj{X) to Y, sending (xq, . . . , Xp) 
to the common value /{xq) = ■ ■ ■ — f(xp), and abusing notation a little we will 
denote this map by / as well. 



5.8.1. The Descent Double Complex. Let C{W^{X)) denote the singular co-chain 
complex of W^{X) (refer to [61] for definition). For each p > 0, we now define a 
homomorphism , 



(5P : G'{W]{X)) ^ C'{WP+\X)) 
as follows: for each i,0 < i < p, define Hp^i : W^{X) W^^'^{X) by 



{'Kp^i drops the «-th coordinate). 

We will denote by (Trp,^)* the induced map on C.(W'}'(X)) ^ C,{W^f^{X)) and 
let TT*p^i : C'{WPf^{X)) ^ G'{W^f{X)) denote the dual map. For e C*{W^j{X)), 
we define 5^ cj) by 



(5.27) <5^0^|](-l)V;+i,, 

1=0 



Definition 5.52 (Descent Double Complex). Now, let D'''{X) denote the double 
complex defined by 'D^'''{X) — {W^ {X)) with vertical and horizontal homomor- 

phisms given by c?* = { — l)^d'^ and 5 respectively, where d is the singular coboundary 



operator, and S is the map defined in (5.271. Also, let 'D^''^{X) = if p < or 
q<0. 
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-> C^(wj(X)) )■ C^(W}(X)) )■ C^(W]{X)) 



> C^(W°j(X)) 



-> C^(Wl(X)) ^ C^(Wj(X)) 



-> C^{W"^{X)) )■ C^{W}{X)) )■ C^{W](X)) 



-> c^(w"i{x)) > d\w}(x)) )■ c°(w'j(,x)) 



Theorem 5.53. |411l24j For any continuous semi- algebraic surjection f : X Y , 
where X and Y are open semi-algebraic subsets o/ R" and R™ respectively (or, 
more generally, for any locally split continuous surjection f ), the spectral sequence 
associated to the double complex D*'*{X) with Ei = W^(D*''{X)) converges to 
H*(C*(F)) =S H*(r). In particular, 

(1) E"^^ =W{W}{X)), and 

(2) i?oo-H*(Tof(D'''(X))) = H*(y). 

We will also need a well-known construction in homotopy theory called homotopy 
colimits, which we define now. 

5.9. Homotopy Colimits. Let A — {Ai, . . . , An}, where each Ai is a sub-complex 
of a finite CW-complex. Let A[„] denote the standard simplex of dimension n — 1 
with vertices in [n]. For / C [n], we denote by A/ the (^/ — l)-dimensional face 
of A[„] corresponding to /, and by Aj (resp. A^) the CW-complex Ai (resp. 

(J Ai). The homotopy colimit, hocolim(^), is a CW-complex defined as follows. 

Definition 5.54 (Homotopy colimit). 

(5.28) hocolim(yl) = |J A/ x - 

IC[n] 

where the equivalence relation ~ is defined as follows. For / C J C [n], let sj j : 
A/ ^ Aj denote the inclusion map of the face A/ in Aj, and let ijj : Aj ^ Aj 
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denote the inclusion map of Aj in Aj. Given (s, x) E Aj x Aj and {t, y) e Aj x Aj 
with I <Z J , then (s,x) '--^ if and only if t = sij{s) and x = iij{y). 

We have an obvious map 

(5.29) fA ■■ hocohm(yl) — > colim(yl) = ^1"! 

sending (s,x) i-^ a;. Notice that for each x G A^"!, 



/^^(x) = |A,J, 

where /a: = {i | x S Ai\. In particular, f^^{x) is contractible for each x G f^^{x), 
and it follows from the Smale-Vietoris theorem [63 that 

Lemma 5.55. The map /_4 is a homotopy equivalence. 

For > 0, we will denote by hocolim<£(^) the subcomplex of hocolim(^) defined 

by 

(5.30) hocolim<f(y^) = |J A/ x ^7/ ~ 

/C[n],#/<f+2 

The following theorem is the key ingredient in the algorithm for computing Betti 
numbers of arrangements described in Section |8] 

Theorem 5.56. For < j < £ we have, 

HJ(hocolim<^(y^)) ^ (hocolim(y^)) = ^(^["1). 

Proof. By Lemma |5.55| we have that 

(5.31) ff(hocolim(^)) ^ff (^["l), i > 0. 

We also have by construction that the {£ + l)-st skeletons of hocolim<^(^) and 
hocolim(^) coincide, which implies that 

(5.32) (hocolim<£(y^)) = HJ(hocolim(yl)), 0<j<t 

The theorem now follows from ( 5.31[ ) and (5.321 above. □ 



6. Algorithms for Computing the First Few Betti Numbers 

We are now in a position to describe some of the new ideas that make possible the 
design of algorithms with single exponential complexity for computing the higher 
Betti numbers of semi-algebraic sets. 

6.1. Computing Covers by Contractible Sets. One important idea in the al- 
gorithm for computing the first Betti number of semi-algebraic sets, is the con- 
struction of certain semi-algebraic sets called parametrized paths. Under a certain 
hypothesis, these sets are semi-algebraically contractible. Moreover, there exists an 
algorithm for computing a covering of a given basic semi-algebraic set, S C R'^, by 
a single exponential number of parametrized paths. 



ALGORITHMIC SEMI-ALGEBRAIC GEOMETRY AND TOPOLOGY 



43 




6.1.1. Parametrized Paths. We are given a polynomial Q € R[Xi, . . . , Xk] such that 
2{Q,R'') is bounded and a finite set of polynomials V C D[Xi, . . . ,Xk]- 

The main technical construction underlying the algorithm for computing the first 
Betti number in [20J, is to to obtain a covering of a given T'-closed semi-algebraic 
set contained in Z((5,R'^) by a family of semi-algebraically contractible subsets. 
This construction is based on a parametrized version of the connecting algorithm: 
we compute a family of polynomials such that for each realizable sign condition 
a on this family, the description of the connecting paths of different points in the 
realization, 7?.(ct, Z(Q, R*"')), are uniform. We first define parametrized paths. A 
parametrized path is a semi-algebraic set which is a union of semi-algebraic paths 



having the divergence property (see Section 4.3.11. 
More precisely. 

Definition 6.1 (Parametrized paths). A parametrized path 7 is a continuous semi- 
algebraic mapping from V C R*^^^ R*^, such that, denoting hy U — TTi,,,k(y) C 
R*^, there exists a semi-algebraic continuous function £ : U —> [0,+oo), and there 
exists a point a in R*^, such that 

(1) V = {{x,t) \ xeU,0<t< i{x)}, 

(2) yxeU, 7(2;, 0) = a, 

(3) Vx e U, -/lx,i{x)) = X, 
(4) 

Vxeuyye uys e [o,e{x)],yt£ [o,£{y)] 

(7(2;, s) = 7(y,i) ^ s = t), 

(5) 

Vx e C/,V y e C/,V s e [0,min(^(a;),^(?;))] 
(7(x, s) = 7(2/, s) ^ V t < s 7(0;, t) = 7(2/, t)) . 
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Given a parametrized path, 7 : y ^ R , we will refer to U = Tri...k{V) as its 
base. Also, any semi-algebraic subset U' C U oi the base of such a parametrized 
path, defines in a natural way the restriction of 7 to the base U' , which is another 
parametrized path, obtained by restricting 7 to the set V C V, defined by V = 
{{x,t) \ xeU',0<t< £{x)}. 

The following proposition which appears in I20j describes a crucial property of 
parametrized paths, which makes them useful in algorithms for computing Betti 
numbers of semi-algebraic sets. 

Proposition 6.2. [20 Let ^ : V ^ be a parametrized path such that U ~ 
7ri...fe(V^) is closed and bounded. Then, the image of j is semi- algebraically con- 
tractible. 

It is also shown in [50] that, 

Theorem 6.3. Moreover, there exists an algorithm that takes as input a finite set 
of polynomials V C R[^i, . . • , X^], and produces as output, 

• a finite set of polynomials A C R[A^i, . • . 

• a finite set O of quantifier free formulas, with atoms of the form P — 0, P > 
0, P < 0, P G A, such that for every semi- algebraically connected compo- 
nent S of the realization of every weak sign condition on V on Z((5,R'^), 
there exists a subset 0(5) C O such that S — |^ 7?.(0, Z((5, R'')), 

6»ee(S) 

• for every Cz Q, a parametrized path 

7e : -> R^ 

with base Ug — 7?.(6', Z(Q, R*^)), such that for each y G 7?,(6', Z(Q, R*^)), 
Im 7e(y, ■) is a semi- algebraic path which connects the point y to a distin- 
guished point ag of some roadmap RM(Z(7'' U {Q},R'')) where V C V, 
staying inside TZ{a{y),Z{Q,K")). 
Moreover, the complexity of the algorithm is s'' +ic?o(fc )^ inhere s is a bound on 

the number of elements ofV and d is a bound on the degrees of Q and the elements 

ofV. 

6.1.2. Constructing Coverings of Closed Semi-algebraic Sets by Closed Contractible 
Sets. The parametrized paths obtained in Theorem |6.3| are not necessarily closed 
or even contractible, but become so after making appropriate modifications. At the 
same time it is possible to maintain the covering property, namely for any given 
T'-closed semi-algebraic 5 set, there exists a set of modified parametrized paths, 
whose union is 5. Moreover, these modified sets are closed and contractible. We 
omit the details of this (technical) construction referring the reader to [5D] for more 
detail. Putting together the constructions outlined above we have: 

Theorem 6.4. There exists an algorithm that given as input a V-closed and 
bounded semi- algebraic set 5, outputs a set of formulas {4>i, . . . , iPm} such that 

• each TZ{(l>i,Yy^) is semi- algebraically contractible, and 

• U 7^(0„R'^■) =Ext(5,R'), 

l<j<M 

where R' is some real closed extension of R. The complexity of the algorithm is 
bounded by s'^^^^^ d^^^ \ where s = fj^V and d = maxpg-p deg(P). 
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6.2. Computing the First Betti Number. It is now an easy consequence of 



the existence of single exponential time covering algorithm (Theorem 6.4 1, and 
Theorem |5.36| stated above, along with the fact that we can compute descriptions 
of the connected components of semi-algebraic sets in single exponential time, that 
we can compute the first Betti number of closed and bounded semi-algebraic sets in 



single exponential time (see Remark 5.37 above), since the dimensions of the images 
and kernels of the homomorphisms of the complex, L2(C(>5')) in Theorem 5.36 can 
then be computed using traditional algorithms from linear algebra. As mentioned 
earlier, for arbitrary semi-algebraic sets (not necessarily closed and bounded), there 
is a single exponential time reduction to the closed and bounded case using the 
construction of Gabrielov and Vorobjov [20l HQ] . 

6.3. Computing the Higher Betti Numbers. 

6.3.1. Double Complexes Associated to Certain Covers. We now describe how cov- 
ers by contractible sets of a closed and bounded semi-algebraic set, S, can be used 



for computing the higher (than the first) Betti numbers of S. Recall (Remark 5.39) 
that it is no longer possible to use the nerve complex of a (possibly non-Leray) 
cover by contractible sets for this purpose. 

In this section, we consider a fixed family of polynomials, V C R[Xi, . . . ,Xk], 
as well as a fixed P-closed and bounded semi-algebraic set, S" C R'^. We also fix a 
number, £,0 < £ < k. 

We define below (see Definition |6.5| below) a finite set of indices. As, which we 
call the set of admissible indices, and a map that associates to each a G As a 
closed and bounded semi-algebraic subset Xa C S, which we call an admissible 
subset. The reason behind having the set of indices is that in the construction of 
our complex the same set might occur with different indices and we would like to 
distinguish these occurrences from each other. 

To each a S Ag, we will associate its level, denoted level(a), which is an integer 
between and £. The set A5 will be partially ordered, and we denote by an(a) C 
A5, the set of ancestors of a under this partial order. For a, (3 G Ag, (3 G an(a) 
will imply that C Xp. 

For each admissible index a £ As, we define a double complex, M'''{a), such 
that 

(6.1) ff(Tof(M*'*(a))) =^ ff(X„), 0<i<£- level(a), 
and for each pair a,/3 G As with a G an(/3) a homomorphism, 

(6.2) r'^'p:M'''ia)^M'-'(,P), 

which induces the restriction homomorphisms between the cohomology groups via 
the isomorphisms in 

<^:ff(X„)^ff(X^) 



for < i < £ — level(a) via the isomorphisms in (6.1 1. 

The main idea behind the construction of the double complex A4*'*{a) is a re- 
cursive one. Associated to any cover of Xa there exists a double complex (the 
Mayer- Vietoris double complex) arising from the generalized Mayer- Vietoris ex- 



act sequence (see Section 5.6 1. If the individual sets of the cover of X are all 



contractible, then the first column of the 'i?i-term of the corresponding spectral 



sequence (cf. Eqn. (5.23 1) is zero except at the first row. The cohomology groups 
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of the associated total complex of the Mayer- Vietoris double complex are isomor- 
phic to those of and thus in order to compute 5o(Xa), . . . , ^f_level(Q)("'^")' 
suffices to compute a suitable truncation of the Mayer-Vietoris double complex. 
Computing (even the truncated) Mayer-Vietoris double complex directly within a 
single exponential time complexity is not possible by any known method, since 
we are unable to compute triangulations of semi-algebraic sets in single exponen- 
tial time. However, making use of the cover construction recursively, we are able 
to compute another double complex, M*'*{a), which has much smaller size, but 
whose associated spectral sequence, 'i?*, is isomorphic to the one corresponding to 



the Mayer-Vietoris double complex. Hence, by Theorem 5.43 (Comparison The 



orem) Tot*(A^*'*(Q;)), is quasi-isomorphic to the associated total complex of the 



Mayer-Vietoris double complex (see Proposition 6.8 below). The construction of 
M.*'*{a) is possible in single exponential time since the covers can be computed in 
single exponential time. 

Finally, given any closed and bounded semi-algebraic set X C R*^ , we will denote 
by C'{X), a fixed cover of X (we will assume that the construction implicit in 



Theorem 6.4 provides such a cover). 

We now define A5, and for each a E A5 a cover C{a) of Xa obtained by enlarging 
the cover C'{Xa). 

Definition 6.5 (Admissible indices and covers). is defined by induction on 
level. 

(1) Firstly, e As, level(O) = 0, Xq = S, an(0) = 0, and C(0) = C'(S'). 

(2) The admissible indices at level i + I are now inductively defined in terms 
of the admissible indices at level < i. 

The set of admissible indices at level i 4- 1 is 

(6.3) IJ {ao-ai---aj\aieC{a),0<j<l-i + l}, 

where IJ denotes the disjoint union, and for each P = aQ ■ ai ■ ■ ■ aj we set 

^13 = ^ao n • • • n Xaj ■ 

We now enlarge the set of ancestor relations by adding: 

(a) For each {ao, . . . ,am} C {/3o, ■ • ■ ,/3n} C C{a), with n < £ — i + 1, 
ao ■ • • Q!m e an(/3o ■ ■ ■ /3„), and a e an(/3o ■ ■ ■ /3„). 

(b) Moreover, if ai ■ ■ ■ am-, Po ' ' ' f^n G A5 are such that for every j € 
{0, . . . , n} there exists i G {0, . . . , m} such that is an ancestor of (3j, 
then ao • • • am is an ancestor of /3o • • • f3n ■ 

(c) The ancestor relation is transitively closed, so that ancestor of an 
ancestor is also an ancestor. 

Finally, for each a E As at level i -I- 1, we define C{a) as follows. Let 
an(a) = {ai, . . . , aAr}. 

(6.4) C{a)=\JC'{Pi---(3N-a). 

where the disjoint union is taken over all tuples (/3i, . . . , Pn) satisfying for 
each 1 < i, j < N , Pi E C{ai), and if E an(Q;j) then Pi E an(/3j). 

Observe that by the above definition, if a,/? G A5 and P E an(a), then each 
a' E C{a) has a unique ancestor in each C(/3), which we will denote by aa,f3{a'). 
The mappings aa^p has the property that ii P E an(a) and 7 E an(/3), then 

aQ,7 — a^,7 ° O'a.,13- 
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Now, suppose that there is a procedure for computing C'(X), for any given V'- 
closed and bounded semi-algebraic set, X, such that the number and the degrees 
of the polynomials appearing the descriptions of the semi-algebraic sets, Xa,a G 
C'{X), is bounded by 

(6.5) D*^", 

where ci > is some absolute constant, and D — J^peV '^'^si^)- 

Using the above procedure for computing C'{X), and the definition of A5, we 
have the following quantitative bounds on and the semi-algebraic sets Xa , a € 
As, which is crucial in proving the single exponential complexity bound of the 
algorithm for computing the first few Betti numbers of semi-algebraic sets. 

Proposition 6.6. [TO] Lei S" C R*" be a V -closed semi- algebraic set, where V C 
R[Xi, . . . , Xk] is a family of s polynomials of degree at most d. Then ^As, as well 
as the number of polynomials used to define the semi-algebraic sets Xa, a G A5 and 
the the degrees of these polynomials, are all bounded by (sd)'^ * 

6.3.2. Double Complex Associated to a Cover. Given the different covers described 
above, we now associate to each a G As a double complex, A4*'*{a), and for every 
(3 € As, such that a G an(/3), and level(a) = level(/3), a restriction homomorphism: 

r'a',:M'''{a)^M'-'{P), 

satisfying the following: 
(1) 

(6.6) H'(Tof(Al*'*(a))) ^ (^a), for < i < ^ - level(a). 



(2) The restriction homomorphism 



r 



M'''{a) M*'*{f3), 



induces the restriction homomorphisms between the cohomology groups: 
r:^^:ff(X„)^ff(X^) 



for < i < £ — level(a) via the isomorphisms in (6.6 1. 

We now define the double complex Ai''*{a). The double complex A4*''{a) is 
constructed inductively using induction on level(a). 

Definition 6.7. The base case is when level(a) = £. In this case the double 
complex, M*'*{a) is defined by: 

M"^'{a) = e„„,cM H°(X„J, 
MP''^{a) = 0, if(7>0orp> 1. 



This is shown diagrammatically below. 
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-> 



The only non-trivial homomorphism in the above complex, 

a;oGC(a) ao>«iGC(a) 



is defined as follows. 

For every P e As, such that a e an(^), and level(a) = level(/3) = we define 
^ajJ M^'^ia) ^ M^'°{(3), as follows. 

Recall that, A^°'° (a) = H°(X„o), and = H°(X;3j. 

ao e C(a) /3o 6 C(/3) 

For (/) e A^°'°(a) and /3o e we define 



We define r^'° : ^ A^^'°(/3), in a similar manner. More precisely, for 

(t> e M°'°{a) and /3o,/3i e C(/3), we define 



(The inductive step) In general the A^P''(q;) are defined as follows using induction 
on level(a) and with Ua = i — level(a) + 1. 



7W0'«(a) ==0, < g, 

MP^Ha) = e„„<...<„^, „,ec(a) Tot«(M*'*(ao • • • a^)), < p, < p + g < 

7WP'«(a) = 0, else. 



The double complex A4*'*{a) is shown in the following diagram: 
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T 

_^ Tot"°-l{M*'*(ao •«!)) 



0«»1 

d-t d 

^ © Tot"«-2(M*'*(ao-ai)) ... q 

<»0<<*1 



^ © Tot2(M'.'(ao -ai)) ... Q 

^ © Totl(Al''*(ao-ai)) ... ^ 

4 ° ' d^- 4 

© hO(S„o) £ Tot°(M''*(ao-ai)) ... © Tot" (M*'*(ao ■•■«„„ )) 

c.oeC(X) ao«»l oo<- ■<»n„ 

The vertical homomorphisms, d, in A^*'* (a) are those induced by the differentials 
in the various 

Tot'{M'''{ao ■ ■ ■ ap)),a^ S C(a). 
The horizontal ones are defined by generalized restriction as follows. Let 

cj>e Tot«(Af'*(ao---ap)), 

ao<---<ap,aiEC{a) 

with 

0<J<9 

and 

'/'io,...,a,eM^''-^"(ao---ap). 

We define 



S: Tot«(A^*'*(ao---ap)) ^ Tot\M'''{ao ■ ■ ■ a^+i)) 

ao<---<ap,aieC{a) ao<---<ap+i 

by 

'5(0)a;ov,ap+i = ( — 1) ''ao---tfr--ap+i,ao---ap+i (^ao,--,c«i,.-,ap+i)' 

0<i<p+l 0<j<q 

noting that for each i,0 < i < p + I, ao ■ ■ • di ■ ■ • a^+i is an ancestor of ao • • • ap+i, 
and 

level(ao • • • dj • • • ctp+i) = 1cvc1(q:o • • • ctp+i) = 1cvc1(q!) + 1, 
and hence the homomorphisms ri''^ i- „ „ „ arc already defined by induc- 
tion. 

Now let, a,(3 £ As with a an ancestor of /? and level(a) = level(/?). We define 
the restriction homomorphism, 

as follows. 
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As before, for (j) e M°'°{a) and /3o G C(/3) we define 
For 0<p, 0<p+g<^ — level(Q;) + 1, we define 
as folfows. 

Let e 7WP'«(a) = Tot«(7W^*(ao ••• ap)). We define 

Qo<---<Qp. ai£C(a) 

'^a'Ic.(/3o---/3p),/3o---/3p'^a^.„(/3o),...,a^,„(/3p); 

/3o<---</3p./3iGC(/3)0<i<9 

where a/3,Q(/3o • ■ ■ /3p) = ap,a{PQ) ■ ■ ■ ap^a{Pp)- Moreover, 

level(a;3,a(/3o • • • /3p)) = level(/3o ■ ■ ■ j3p) = level(Q:) + 1, 
and hence we can assume that the honiomorphisms r*'* „ ^ a a used in the 
definition of r*'^ are aheady defined by induction. 

The foUowing proposition proved in [10 is the key ingredients in the single 
exponential time algorithm for computing the first few Betti numbers of semi- 
algebraic sets. 

Proposition 6.8. [lOj For each a G A5 the double complex A^*'*(a) satisfies the 
following properties: 

(1) (Tot* (>!•'• (a))) ^ WiXa) forO<i<£~ level(a). 

(2) For every (3 S A5, such that a is an ancestor of (3, and level(Q;) — level(/?), 
the homomorphism, r*'* : A^*'*(a) ^A'''{P), induces the restriction 
homomorphisms between the cohomology groups: 

for < i < £ — level(Q!) via the isomorphisms in (1). 
Proof Sketch. For the benefit of the reader we include an outline of the proof of 



Proposition 6.8 referring the reader to |10j for more detail. The proof is by induction 
on level(a) . After having chosen a suitably fine triangulation A of 5 which respects 
all the admissible subsets Xq., we construct by induction on level(Q;), for each 
a £ A5, a double complex D*'*(a) and homomorphisms 

(j)''' : A<*'*(a) ^D*'*(a), 

and 

V-* :C*(A„)^Tof(D*'*(a)), 

where Aq, denotes the restriction of the triangulation A to X^- It is then shown (in- 
ductively) that each homomorphism in the following diagram is a quasi-isomorphism. 

Tof(D*'*(a)) 

(6.7) 




Tof(7W'*(a)) 

This suffices to prove the proposition. □ 
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We now give an example of the construction of the complex described above in 
a very simple situation. 

Example 6.9. We take for the set S, the unit sphere C R'^. Even though 
this example looks very simple, it is actually illustrative of the main topological 
ideas behind the construction of the complex A4*'*(S) starting from a cover of S 
by two closed hemispheres meeting at the equator. Since the intersection of the 
two hemisphere is a topological circle which is not contractible, Theorem 5.34 is 
not applicable. Using Theorem 5.36 we can compute H"(S'), H^(S'), but it is not 
enough to compute H^(S'). The recursive construction of A4''* described in the 
last section overcomes this problem and this is illustrated in the example. 




Figure 10. Example of C R'^ 



We first fix some notation (see Figure 10 1. Let Hi and H2 denote the closed 
upper and lower hemispheres respectively. Let H12 = Hi D H2 denote the equator, 
and let H12 — Ci U C2, where Ci,C2 are closed semi-circular arcs. Finally, let 
C12 = Ci n C2 — {Pi,P2}, where Pi,P2 are two antipodal points. 

For the purpose of this example, we will take for the covers C the obvious ones, 
namely: 



C'{S) 


= {HuH2}, 




cm 


= 


z = l,2, 


C'{Hi2) 


= {CuC2}, 






= m, 


z = 1,2, 


C'{Ci2) 


= {PuP2}, 




C'iPr) 


= {P^}, 


i = 1,2. 



Note that, in order not to complicate notation further, we are using the same 
names for the elements of C'(-), as well as their associated sets. Strictly speaking, 
we should have defined, 

C'{S) = {ai, 02}, Xa^ ^ Hi,Xa2 = H2, 

However, since each set occurs at most once, this does not create confusion in this 
example. 
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Note that the elements of the sets occurring on the right are all closed, bounded 
contractible subsets of S. It is now easy to check from Definition |6.5| that the 
elements of As in order of their levels as follows. 



(1) Level 0: 



and 



e As,level(0) = 0, 



C(0) = {ai,a2}, Xai^Hi,Xa2 

(2) Level 1: The elements of level 1 are 

ai,a2,ai ■ a2, 

and 

C{a2) - {/32}, Xp,=H2, 
C{ai-a2) = {03,^4}, Xp.^=Ci,Xp^^C2. 

(3) Level 2: The elements of level 2 are f3i, (32, f33, P4, - Pa- We also have, 



cm 

CiPs • Pi) 



{7.}, 
{7.}, 

{75,76}, X^ 



X, 



1,2, 



Ci-2, i = 3,4, 
Pi-i, i = 5,6. 



We now display diagrammatically the various complexes, A4*'*{a) for a e Ag 
starting at level 2. 

(1) Level 2: For 1 < i < 4, we have 



> 

H°(X^.) )• 

Notice that for 1 < i < 4, 

HO(Tof(Al-'*(A)))-HO(Xft) = 
The complex Ai''*{f33 ■ fS^) is shown below. 

> 



h°(Fi)©h"(P2) s> 

Notice that, 

HO(Tof(>C^'(/33 • /34))) = H0(X^3.;5j - 
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(2) Level 1: For i = 1, 2, the complex M'''{ai) is as follows. 







-4- 







-> 



Notice that for i = 1,2 and j = 0,1, 

H^(Tof(7W'*(a,))) 
The complex M.*'*{ai ■ a2)is shown below. 











h''(Ci)0H°(C2) > H°(Pi)0H°(P2) 

Notice that for j = 0, 1, 

(Tot* (ai • as))) ^ H^(fl"i2). 

(3) Level 0: 

The complex A^*'*(0) is shown below: 











h°(Pi)0hO(P2) 



H°(Hi)©H°(/f2) > H°(Ci)0H°(C2) 











-4- 




The matrices for the homomorphisms, 6^'^ and in the obvious bases 
are both equal to 

■ 1 1 
1 1 
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From the fact that the rank of the above matrix is 1 , it is not too difficult 
to deduce that, ff (Tof(A^*^*(0))) ^ W{S), for j = 0, 1,2, that is 

HO(Tof(X*^*(0))) ^ Q, 
Hi(Tof(Al*^*(0))) ^ 0, 
H2(Tot*(Al*^*(0))) ^ Q. 

6.3.3. Algorithm for Computing the First Few Betti Numbers. Using the construc- 



tion of the double complex outline in Section 6.3.2 as well as the single exponential 
time algorithm for obtaining covers by contractible sets described in Section 6.1.2 
along with straightforward algorithms from linear algebra, it is now easy to obtain 
the following result: 

Theorem 6.10. |10j For any given I, there is an algorithm that takes as input a V- 
formula describing a semi-algebraic set S C R*^, and outputs bQ(S), . . . , bi{S). The 
complexity of the algorithm is (sd)'^'^* , where s = #(7^) and d = maxpg-p deg(P). 

Note that the complexity is single exponential in k for every fixed £. 

7. The Quadratic Case 

7.1. Brief Outline. We denote by S'^ C R*^^^ the unit sphere centered at the 
origin. Consider the case of semi-algebraic subsets of the unit sphere, C R'^'^"'^, 
defined by homogeneous quadratic inequalities. There is a straightforward reduc- 
tion of the general problem to this special case. 

Let S' C S'" be the set defined on S'" by s inequalities, Pi < 0,...,Ps < 0, 
where Pi, . . . ,Ps G R[-'^'oj ■ ■ ■ ,Xk] are homogeneous quadratic polynomials. For 
each i,l < i < s, let Si C S*" denote the set defined on S'' by Pi < 0. Then, 

s 

S^f]S,. There are two main ingredients in the polynomial time algorithm for 

i=l 

computing the top Betti numbers of S. 

The first main idea is to consider S as the intersection of the various Si's and 
to utilize the double complex arising from the generalized Mayer- Vietoris exact 
sequence (see Section [s]) . It follows from the exactness of the generalized Mayer- 



Vietoris sequence (see Definition 5.44), that the top dimensional cohomology groups 
of S are isomorphic to those of the total complex associated to a suitable truncation 
of the Mayer- Vietoris double complex. However, computing even the truncation 
of the Mayer- Vietoris double complex, starting from a triangulation of S would 
entail a doubly exponential complexity. However, we utilize the fact that terms 
appearing in the truncated comple x dep end on the unions of the SiS taken at 

most £ + 2 at a, time (cf. Remark 



5.501. There are at most ( . j such sets 



Moreover, for semi-algebraic sets defined by the disjunction of a small number of 
quadratic inequalities, we are able to compute in polynomial (in k) time a complex, 
whose homology groups are isomorphic to those of the given sets. The construction 
of these complexes in polynomial time is the second important ingredient in our 



algorithm and is outlined below in Section 7.4 These complexes along with the 
homomorphisms between them define another double complex whose associated 
spectral sequence (corresponding to the column-wise filtration) is isomorphic from 
the E2 term onwards to the corresponding one of the (truncated) Mayer- Vietoris 
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double complex (see Theorem 7.21| below). Since, we know that the latter converges 



to the homology groups of S, the top Betti numbers of S are equal to the ranks 
of the homology groups of the associated total complex of the double complex we 
computed. These can then be computed using well known efficient algorithms from 
linear algebra. 

In order to carry through the program described above we need to understand 
a few things about the topology of sets defined by homogeneous quadratic inequal- 
ities. 

7.2. Topology of Sets Defined by Quadratic Inequalities. In this section we 
state a few results concerning the topology of sets defined by quadratic inequali- 
ties, which are exploited in designing efficient algorithms for computing their Betti 
numbers. 

7.2.1. Case of One Quadratic Form. We first consider the case of a single quadratic 
form Q G R[Xo, . . . ,Xk]- Let S* C S*^ be the set defined by Q > on the unit sphere 
in R*^^^. The crucial fact that distinguishes quadratic forms from forms of higher 
degree is that the homotopy type of the set S is determined by a single invariant 
attached to the quadratic form Q, namely its index. 

Definition 7.1 (Index of a quadratic form). For any quadratic form Q, mdex{Q) 
is the number of negative eigenvalues of the symmetric matrix of the corresponding 
bilinear form, that is of the matrix M such that, Q{x) = {Mx,x) for all x G R*^^^ 
(here (•, •) denotes the usual inner product). We will also denote by Xi{Q),0 < i < k, 
the eigenvalues of M , in non-decreasing order, i.e. 

Ao(Q) < Ai(Q) < ••• < Afe(Q). 

A simple argument involving diagonalizing the quadratic form Q (see [5] or [TT]) 
yields that the homotopy type of the set S defined above is related to the index(Q) 
by 

Proposition 7.2. The set S is homotopy equivalent to the s''^'"'*°''('3) . 



Example 7.3. The following figure (Figure 111 illustrates Proposition 7.2 We 
display (from left to right) the subsets of described by the inequalities 

+ 2X1 + 3^1 > (index = 0), 

Xl + 2X1 - 3^1 > (index =1), 

Xl - 2X1 - 3^2 > (index = 2) 

respectively. Notice that each of the quadratic forms defining these sets are already 
in a diagonal form, and hence its index can be read off directly from the signs of 
the coefficients (the index is the number of negative coefficients) . By Proposition 
7.2 these sets have the homotopy types of S^, S^, and S° respectively, as can be 



also seen from the displayed images below. 





Figure 1 1 . Subsets of defined by one homogeneous quadratic 
inequality of index 0, 1 and 2. 

7.2.2. Case of Several Quadratic Forms. Now let Qi, . . . , Qs be homogeneous qua- 
dratic polynomials in R[Xo, . . . , Xi^]. 

We denote by Q = (Qi, ...,Qs) : R^^^ ^ R", the map defined by the forms 
Qi, - ■ ■ , Qs- 

Let 

(7.1) T= \J {xeS'' \Q,ix) <0}, 

l<i<s 

and let 

(7.2) n^lLueR" \ \uj\^ l,uj,< 0,1 <i<s}. 
For cij G we denote by ujQ the quadratic form 

s 

(7.3) wQ = y^uj^Qj. 

1=1 

Let S C X S'' be the set defined by 

(7.4) B = {{uj, x)\uj en,x and ujQ{x) > 0}. 

We denote hy 4>i : B ^ and (j)2 B ^ S'^ the two projection maps. 

B 

4>i / \ 4)2 



With the notation developed above we have 

Proposition 7.4. The map 02 gives a homotopy equivalence between B and 
MB) = T. 

We denote by 

(7.5) nj = {Luen I Xj{ujP)>0}. 

It is clear that the fij's induce a filtration of the space H., i.e., fio C fii C • • • C 

^k- 

The following lemma follows directly from Proposition |7.2[ It is an important 
ingredient in the algorithms for computing the Euler-Poincare characteristic as 
well as the Betti numbers of semi-algebraic sets defined by quadratic inequalities 
described later. 

Lemma 7.5. The fiber of the map (pi over a point uj £ ilj-i has the homotopy 
type of a sphere of dimension k — j. 

We illustrate Lemma |7.5| with an example. 
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Example 7.6. In this example s = 2,k — 2, and 



]i ^ - X„ - - X. 



2 7 

2 



^ ^ o o o o 



Figure 12. Type change: 



is not shown. 



The set is the part of the unit circle in the third quadrant of the plane. In 
the following Figure 12 we display the fibers of the map (p'^^{Lij) C i? for a se- 
quence of values of w starting from (—1,0) and ending at (0,-1). We also show 
the spheres of dimensions 0, 1, and 2, that these fibers retract to. At u; — (—1,0), 
it is easy to verify that index{ujQ) — 3, and the fiber (p^^{uj) C B is empty. 
Starting from uj = (— cos(arctan(l)), — sin(arctan(l))) we have index(a;(5) = 2 
and the fiber Lp^^(iS) consists of the union of two spherical caps homotopy equiv- 
alent to S°. Starting from cj = (— cos(arctan(l/2)), — sin(arctan(l/2))) we have 
index(a;(5) — 1, and the fiber (^^^(cj) is homotopy equivalent to S^. Finally, start- 
ing from Lo — (— cos(arctan(l/3)), — sin(arctan(l/3))), index(a;Q) = 0, and the 
fiber Lp^^iuS) stays equal to to S^. 



As a consequence of Lemma 7.5 we obtain the following proposition which relates 



the Euler-Poincare characteristic of the set T (cf. Eqn. (7.1 1) with the Borel-Moore 
Euler-Poincare characteristic (cf. Definition 5.221 of fij \ rij-i, < j < fc + 1. 



Proposition 7.7. 

fc+i 

(7.6) x(r) - x^'^'^T) ^ V x^*^(% \ f^,-i)(i + (-i)('=-^")). 



J=0 



It is instructive to continue Example |7.6| and compute the Euler-Poincare char- 
acteristic of the set T in that case using Proposition |7.7| 



Example 7.8. In this example for each < < 3, 

\ rij_i is homeomorphic to [0,1) 

and for j = we have 

r^o \ i^-i i^o is homeomorphic to [0, 1]. 
Recall from Eqn. that x^*^([0,l)) = and X^^^([0,1]) 



Finally using Eqn. (7.6) we deduce that 



x([o,i]) 



X{T) - 2. 
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7.3. Computing the Euler-Poincare Characteristics of Sets Defined by 
Few Quadratic Inequalities. Proposition |7. 7| reduces the problem of computing 
the Euler-Poincare characteristic of the set T, which is defined by quadratic forms in 
fc-l- 1 variables, to computing the Borel-Moore Euler-Poincare characteristics of the 
sets flj \ ilj-i C il, which are defined by 0{k) polynomials in s variables of degree 
also bounded by 0{k). We now utilize an efficient algorithm for listing the Borel- 
Moore Euler-Poincare characteristics of the realizations of all all sign conditions of 
a family of polynomials developed in |18| to compute the terms occurring on the 
right hand side of Eqn. ( |7.6| . The complexity of this algorithm is exponential in 
the number of variables (which is 0{s) in this case) and polynomial in the number 
and degrees of the input polynomials (which are 0{k) in this case). 

The set T is defined by a disjunction of homogeneous quadratic inequalities. But 
since the Euler-Poincare characteristic satisfies the inclusion-exclusion formula (cf. 



Proposition 5.251, we are able to compute it for sets defined by a conjunction of 
such inequalities within the same asymptotic time bound by making at most 2^ 
calls to the algorithm for disjunctions. 

For inhomogeneous quadratic inequalities there is an easy reduction to the ho- 
mogeneous case (see [3] for detail) . As a result we obtain 

Theorem 7.9. |9] There exists an algorithm which takes as input a closed semi- 
algebraic set 5 C R ' defined by 

Pi <0,...,P, <0,P, eR[Xi,...,Xfc],deg(F,) <2, 

and computes the Euler-Poincare characteristic of S. The complexity of the algo- 
rithm is fc*^*^") . 

Remark 7.10. Very recently [13] the above algorithm has been generalized to the 
following setting. Let 

Q = {Qi,...,gjcR[Xi,...,Xfe,yi,...,y,] 

with degxiQi) < 2,degy(Qi) < d,l < i < s, and 

r CR[Yi,...,Yi] 

with deg(P) < d,P <E V and #P = to. Let S c R''^^^ be a 7^ U Q-closed semi- 
algebraic set. Then, 

Theorem 7.11. |14j There exists an algorithm for computing the Euler-Poincare 
characteristic of S whose complexity is bounded by {kimd)'^^'^'^^'^^^\ 



Notice that Theorem 7.11| is a generalization of Theorem |7.9| in several respects. 



It allows a subset of the variables to occur with degrees bigger than 2 (and the 
complexity of the algorithm is exponential in the number of these variables) and it 
takes as input general V U Q-closed semi-algebraic sets, not just basic closed ones. 

7.4. Computing the Betti Numbers. 

7.4. L The Homogeneous Case. We first consider the homogeneous case. 

Let V — (Pi,...,Ps) C R[Xo, . . . , Xj.] be a s-tuple of quadratic forms (i.e. 
homogeneous quadratic polynomials). For any subset Q C P we denote by Tq C S'^ 
the semi-algebraic set 

Pq = U ^ I ^(^) ^ 0}, 
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and let 

S= f]{xGS'' I P{x)<0}. 
Per 

We denote by C*{H{Tq)) the co-chain complex of a triangulation 'H{Tq) of Tq 
which is to be chosen sufficiently fine. 

We first describe for each subset Q CV with = £ < k a, complex, A4q, and 
natural homomorphisms 

i^Q--C'{H{TQ))-^M'Q 

which induce isomorphisms 

ra : H*{C'{H{Tq))) ^ H*{M'q). 

Moreover, fov B C AcV with i^A = < k, we construct a homomorphism 

of complexes 

such that the following diagram commutes. 



(7.7) 



H*{c'{n{TjO)) H*ic'{n{Ts))) 

In the above diagram (p^ ^ and r* are the induced homomorphisms of (t)A,B and 
the restriction homomorphism r respectively. 

Now consider a fixed subset Q cP, which without loss of generality, we take to 
be {Pi,..., Pi}. Let 

P=(Pi,...,P,):R'=+^^R^ 

denote the corresponding quadratic map. 
Let r2 = and 

CIq = {wgR^ \ \u}\ = l,Wi<0,l<i<e}. 

Let Bq c f^Q X S*" be the set defined by 

Bq = {{u>,x) \ u> G Qq, X eS'^ and u>P{x) > 0}, 

and wo denote by (pi^Q : Bq and (j)2,Q ■ Bq S'^ the two projection maps. 

For each subset Q' C Q we have a natural inclusion Qq/ ^ Qq. 

7.4.2. Index Invariant Triangulations. We now define a certain special kind of semi- 
algebraic triangulation of ilg that will play an important role in our algorithm. 

Definition 7.12 (Index Invariant Triangulation). An index invariant triangulation 
of consists of: 

(1) A semi-algebraic triangulation, 

h : Aq — > Q.Q 

of fig which is compatible with the subsets Qqi for every Q! G Q and such 
that for any simplex a of Ag, index(wPQ) as well as the multiplicities of 
the eigenvalues of ujPq stay invariant as oj varies over /i((t); 
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(2) for every simplex cr of Ag with index(cLi Pg) ~ j for lu G h{a), a uni- 
form description of a family of orthonormal vectors eo((T, w), . . . ,efc((T, w), 
parametrized by (w, x) e h((j) having the property that 

{ej{a,uj),...,ek{(J,i^)} 

is a basis for the linear subspace L~^{i-o) C R'^"''"'^ (which is the orthogonal 
complement to the sum of the eigenspaces corresponding to the first j 
eigenvalues of t^Po)- 

An algorithm to compute index invariant triangulations is described in (see 
also HU) the complexity of this algorithm is bounded by k"^"^"'. The same bound 
holds for the size of the complex Ag as well as the degrees of the polynomials 
occurring in the parametrized representation of the vectors {eo((T, w), . . . , efe(cr, lo)}. 

Now fix an index invariant triangulation h : Aq ^Iq satisfying the complexity 
estimates stated above. 

We now construct a cell complex homotopy equivalent to Bq . It is obtained by 
glueing together certain regular cell complexes, /C((t), where a e Ag. 



Figure 13. The complex Ag. 

Let 1 ^ £o ^ El ^ • • • ^ £^ > be infinitesimals. For t G Ag we denote by 
Dr the subset of f defined by 

Dr — {v £ f I dist(w, 6) > £dim(6i) for all 9 ~i a} 

where dist refers to the ordinary Euclidean distance. Now let a ^ t be two simplices 
of Ag. We denote by D^^t the subset of f defined by 

Da,T ^ {v £ f I dist(z), cr) < edim(cr)j and dist(w,6') > £dim(6() for all & -< f}. 
Note that 

|Ag|= \J D^U y D,^,. 

Also, observe that the various D^-'s and Dcr^s are all homeomorphic to closed 
balls and moreover all non-empty intersections between them also have the same 
property. 
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Dr n 



Figure 14. The corresponding complex C( A g). 



Definition 7.13. The union of the Dr's and D^^s together with the non-empty 



intersections between them form a regular cell complex (cf. Definition 5.4), C(Ag) 
whose underlying topological space is |Aq| (see Figures Il3] and 14 1. 



We now associate to each Da (respectively, -Dcr.r) a regular cell complex, /C(cr), 
(respectively, JC{a, r)) homotopy equivalent to 4>i^{h{Dcr)) (respectively, 4>i^{h{Da.T))- 

For each a € Aq and uj £ h(a) let {eo(cr, w), . . . ,efc((T, w)} be the continuously 
varying orthonormal basis of R*^^^ computed previously. 

The orthonormal basis 

{e(){a,uj), . . . ,efe(cr,w)} 
determines a complete flag of subspaces, T{a,LiS), consisting of 
i^°(a,a;) = 0, 
F'^{a,uj) = span(efe(cr, w)), 
F'^{a,u},x) = span(efc(CT,w),efc-i(a-,tj)), 



(cr,w) ^ K 



fe+i 



Definition 7.14. For < j < fc let c^(ct,(jj) (respectively, cj(f7,tj)) denote the 
(fc — j)-dimensional cell consisting of the intersection of the F^^^^^{a, oj) with the 
unit hemisphere in R*^^^ defined by 

{x e S*-' I {x,ej{(j,uj)) > 0} 

(respectively, {a; G S*^ | {x,ej{a,uj)) < 0} ). 

The regular cell complex /C(ct) (as well as /C((t, r)) is defined as follows. 

For each v G |Ag| and cr G Ag let w(cr) € |cr| denote the point of |cr| closest to v. 

The cells of JC{a) are 

{(x,w) I X e cf{a,Lu),LU e h{c)} 
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where index(a;Pg) < j < k and c G C(Ag) is either I?o- itself or a cell contained in 
the boundary of D^r ■ 

Similarly, the cells of /C((t, r) are 

{{x,Lu) I X e cf{a,h{v{cr))),v^ h-\uj) E c} 

where index(u;Pg) < j < k and c G C(Ag) is either Da-^r itself or a cell contained 
in the boundary of -Do-.t- 

Our next step is to obtain cellular subdivisions of each non-empty intersection 
amongst the spaces associated to the complexes constructed above and thus obtain 
a regular cell complex, IC{Bq), whose associated space, \IC{Bq)\, will be shown to 
be homotopy equivalent to Bq. 

First notice that |/C((t',t')| (respectively, |/C(ct)|) has a non-empty intersection 
with |/C(cr, t)| only if D^'^t' (respectively, D^r') intersects D^^t- 

Let D be some non-empty intersection amongst the Z^o-'s and D^ys, i.e. Z? is a 
cell of C(Ag). Then, D C |r| for a unique simplex t G Ag and 

with ui ^ (72 ^ ■•• ^ (Tp ^ = T and p < i. 

For each i, 1 < j < p+ 1, let {/o((Ti, w), . . . , fk{f7i, v)} denote a orthonormal basis 
of R'^^^ where 

fj{cTi,v) = limej(crj, h(tv{ai) + (1 - t)u((Ti))),0 < j < k, 

and let J-{ai, v) denote the corresponding flag consisting of 
F^'icT.^v) = 0, 
^^((7^,1;) = span(/fc(crj,w)), 
F^{ai,v) = span(/fc(CTi,w),/fe_i(crj,u)), 



F'=+i(f7„i;) = R'=+^ 

We thus have p + 1 different flags 

T{ai,v), . . .,T{crp+i,v), 

and these give rise to p-l- 1 different regular cell decompositions of S*^. 

There is a unique smallest regular cell complex, IC'{D,v), that refines all these 
cell decompositions whose cells are the following. Let L C R*^^^ be a linear sub- 
space of dimension j, < j < fc + 1, which is an intersection of linear subspaces 
Li, . . . , Lp+i where Li G T{(Ji, v),l <i<p+l<£+l. The elements of the flags 
J^((Ti , w) , . . . , J-{ap-i-i ,v) of dimension j + 1 partition L into polyhedral cones of var- 
ious dimensions. The intersections of these con es with S*^ over all such subspaces 

illustrates the refinement described 



15 



L C R''+^ are the cells of K.'{D,v). Figure 

above in case of two fiags in R"^. We denote by /C(D, v) the sub-complex of JC'{D, v) 
consisting of only those cells included in L{ai, h{v{ai))) fl S*". 



We now triangulate h{D) using the algorithm implicit in Theorem 4.5 (Triangu- 
lation) so that the combinatorial type of the arrangement of flags 

T{ai,v), . . .,T{ap+i,v) 
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■?^(<J2,l') 



Figure 15. The cell complex /C'(£',w). 



and hence the cell decomposition IC'{D,v) stays invariant over the image, hoiO), 
of each simplex, 9, of this triangulation. Notice that the combinatorial type of 
the cell decomposition JC'{D,v) is determined by the signs of the inner products 
{fj(.(^i,v),fj'{(T,,,v)) where < < 1 < i, i' < p + 1. 

We compute a family of polynomials Ad C R[Zi,...,Zf] whose signs deter- 
mine the vanishing or non- vanishing of the inner products {fj{(7i, v), fj'{ai', v)),0 < 
j, f ^ k,l < < p+ 1. It is then clear that the combinatorial type of the cell de- 
composition /C'(Z), v) will stay invariant as lu varies over each connected component 
of any realizable sign condition on C R[^i, . . . , Zi]. 

Given the complexity bounds on the rational functions defining the orthonormal 
bases {eo(cr, w), . . . , ee{(J,uj)} lu G h^a), stated above that the number and degrees 
of the polynomials in the family Ad are bounded by A:^ * ' . We then use the 
algorithm implicit in Theorem 4.5 (Triangulation) with Ad as input, to obtain the 
required triangulation. 

The closures of the sets 

{{x,Lu) I xeceic{D,h-^{cj)), uj&h{hD{e))} 

form a regular cell complex which we denote by K,{D). 

The following proposition gives an upper bound on the size of the complex IC{D). 
We use the notation introduced in the previous paragraph. 

Proposition 7.15. For each u G h{D), the number of cells in IC{D,h^^{uj)) is 
bounded by k^^^\ Moreover, the number of cells in the complex IC{D) is bounded 
by fc^°'*' . 

Note that there is a homeomorphism iD,iJi '■ |A^(o'i,T)| n <j)i^{h{D)) — > 
which takes each cell of \IC{ai, t) | n (p^^{h{D)) to a union of cells in IC{D). We use 
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these homeomorphisms to glue the cell complexes IC{ai, t) together to form the cell 
complex ]C{Bq). More precisely 

Definition 7.16. K.{Bq) is the union of all the complexes JC{D) constructed above, 
where we use the maps in.ui to make the obvious identifications. 

We have that 

Proposition 7.17. |/C(i3Q)| is homotopy equivalent to Bq. 
We also have 

Proposition 7.18. The number of cells in the cell complex JC{Bq) is bounded by 

Proposition 7.19. [TP For Q < i < k — 1, the induced homomorphisms 

^^:ff(C-(H(rQ)))^ff(A^-Q) 

are isomorphisms. 

Now XetBcAcV with #^ = #B + 1< fc. 

The simplicial complex Ag is a subcomplex of and hence, Cg* is a subcom- 
plex of C*^* and thus there exists a natural homomorphism (induced by restriction) 

(t^Afi ■ C^' C'b * 

and let 

0^,6 : Tot*(C;^*) = M\^Ml = Tof(C;'*), 

be the induced homomorphism between the corresponding associated total com- 
plexes. 

The complexes A^^,A^g, and the homomorphisms, (/'^,b, i/'^i V'B satisfy 
Proposition 7.20. [TT] The diagram 

M\ — > M'n 



(7.8) i'A 



■<I>B 



C'{n{T^)) C-(H(Tg)) 

is commutative, where r is the restriction homomorphism. 
We denote by 

(l>B,A ■■ Ml M'j^ 

the homomorphism dual to (/)^b. We denote by 2?^'* the double complex defined 
by: 

l^v' = ^ M" 



The vertical differentials 
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are induced component-wise from the differentials of the individual complexes A4q 
The horizontal differentials 

S : VP;" ^ VP+^''' 
are defined as follows: for a E Vlp'^ — (B#q=p+iM'q for each subset 

with io < • • • < ip+i the Q-th component of Sa G D^j^^''^ is given by 
where = Q \ 

— y e#Q=iAl| 
— > (B#q=iMq 

1" 



We have the following theorem. 
Theorem 7.21. ^ For < i < k. 

H'{S) ^H' {Tot' {V'^')). 

Finally, using Theorem |7 . 2 1 1 we have 

Theorem 7.22. [llj There exists an algorithm which given a set of s polynomi- 
als, V = {Pi,...,Ps} C R[Xi, . . . , Xfc], with deg(Pi) < 2, 1 < i < s, computes 
bk-i{S), . . . , bk-e{S), where S is the set defined by Pi < 0, . . . , < 0. The com- 
plexity of the algorithm is 

(7.9) |^M,.o...,., 

i=0 ^ ^ 

// the coefficients of the polynomials in V are integers of bit-sizes bounded by t, 
then the bit-sizes of the integers appearing in the intermediate computations and 
the output are bounded by T^sk)'^"'' * ' 

For certain applications we need the following more detailed version of Theorem 

OH 

Theorem 7.23. [11 There exists an algorithm which takes as input a family of 
polynomials {Pi, . . . jP,} C R[^i . . . with deg(Pi) < 2 and a number £ < k, 



a<j<p+i 



-)#q=2Mq e#Q=3A^S 



j. 1 



d 
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and outputs a complex V''' . The complex Tot*(I?*'*) is quasi-isomorphic to 0^(5*), 
the truncated singular chain complex of S , where 

S= Pi {a; G I P{x) < 0}. 
Pev 

Moreover, given a subset V' <ZV with 

S' = []{x(^ R*-' I P{x) < 0}. 
Pev 

the algorithm outputs both complexes T)*'* and T)'*'* (corresponding to the sets 
S and S' respectively) along with the matrices defining a homomorphism ^-p -pi 
such that : H*(Tof(X>;^*)) ^ H,(5') -> H^S") ^ H*(Tot*(X'';'*)) is the 

homomorphism induced by the inclusion map i : S ^ S' . The complexity of the 
algorithm is X^i^o i'i}^^ ' * 

7.5. Projections of Sets Defined by Quadratic Inequalities. There are two 
main ingredients in the algorithm for computing Betti numbers of projections of sets 
defined by quadratic inequahties. The first is the use of descent spectral sequence 



described in Section 5.8 Notice that the individual terms occurring in the double 
complex in Section |5.8.1| correspond to the chain groups of the fibered products of 
the original set. A crucial observation here is that the fibered product of a set defined 
by few quadratic inequalities is again a set of the same type. However, since there 
is no known algorithm for efficiently triangulating semi-algebraic sets (even those 
defined by few quadratic inequalities) we cannot directly use the spectral sequence 
to actually compute the Betti numbers of the projections. In order to do that we 
need an additional ingredient. This second main ingredient is the polynomial time 



algorithm in Theorem 7.23 for computing a complex whose cohomology groups are 
isomorphic to those of a given semi-algebraic set defined by a constant number 
of quadratic inequalities. Using this algorithm we are able to construct a certain 
double complex, whose associated total complex is quasi-isomorphic to (implying 
having isomorphic homology groups) a suitable truncation of the one obtained from 
the cohomological descent spectral sequence mentioned above. This complex is of 
much smaller size and can be computed in polynomial time and is enough for 
computing the first q Betti numbers of the projection in polynomial time for any 
constant q. 

We have the following theorem. 

Theorem 7.24. [2A\ There exists an algorithm that takes as input a basic semi- 
algebraic set S C K'^"'''" defined by 

Pi>0,...,Pe> 0, 

with Pi e R[Xi, . . . , Xfe, Yi, . . . , Ym], deg(Pi) < 2, 1 < i < £ and outputs 

bo{7T{S)),...,bMS)), 

where tt : M*^^™ M"' be the projection onto the last m coordinates. The complex- 
ity of the algorithm is bounded by {k -\- mf' ' ' . 
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8. Betti Numbers of Arrangements 

In this section we describe an algorithm for computing the Betti numbers of the 
union of a collection, S, of subsets of R*^, where each set is assumed to be a closed 
and bounded semi-algebraic set of constant description complexity. It is customary 
to call the collection S an arrangement and we will refer to this problem as the 
problem of computing the Betti numbers of the arrangement S. A semi-algebraic 
set in R'' is said to have constant description complexity if it can be described by 
a first order formula of size bounded by some constant (see also [12j for a more 
general mathematical framework) . The key point which distinguishes the results in 
this section from those in the previous sections is that unlike before, here we are 
interested only in the combinatorial part of complexity estimates - i.e. the part of 
the complexity that depends on the number of sets in the input. Since the input 
sets are of constant description complexity, the algebraic part of the complexity 
- i.e. the part that depends on the degrees and number of polynomials defining 
each set - is bounded by a constant. This point of view, which is now standard 
in discrete and computational geometry (see (TUSlj), presents new challenges from 
the point of view of designing efficient algorithms for computing Betti numbers of 
arrangements of sets of constant description complexity. 

Notice that, unlike before, in this setting it is not important to obtain a good 
(say single exponential in k) bound on the the algebraic part of the complexity, 
since it is bounded by some constant regardless of the exact nature of the bound. 
Thus, we have much greater flexibility in designing algorithms, since we can utilize 



triangulation algorithms (cf. Theorem 4.5) which have doubly exponential com 



plexity as long as the number of sets in the input to each such call is bounded by 
a constant. 

On the other hand the algorithms described in Section |6] while having single 
exponential complexity, are no longer the best possible in this setting, since the 
combinatorial complexities of these algorithms are very far from being optimal. 
The goal is to use the flexibility afforded in the algebraic part to design algorithm 
having much tighter combinatorial complexity. 

A version of the main result of this section (Algorithm [l] below) appears in [5] 
where a spectral sequence argument is used. We present here a different (and sim- 
pler) algorithm which avoids spectral sequences but instead uses the more geometric 



notion of homotopy colimits (cf. Definition 5.541. The new algorithm has the same 
complexity as the previous one. 

8.1. Computing Betti Numbers via Global Triangulations. As we have seen 
in Section[5] one approach towards computing the Betti numbers of the arrangement 
is to obtain a triangulation of the whole arrangement using the algorithm implicit 
in Theorem |4.5[ Thus, in order to compute the Betti numbers of an arrangement of 



n closed and bounded semi-algebraic sets of constant description complexity in R*^ 
it suffices to first triangulate the arrangement and then compute the Betti numbers 
of the corresponding simplicial complex. However, using the complexity estimate in 
Theorem 



4.5 



the complexity of computing such a triangulation is 0(n^ ). However, 
since the Betti numbers of such an arrangement is bounded by 0(n*) (cf. Theorem 
2.4 1, it is reasonable to ask for an algorithm whose complexity is bounded by 0{n''). 



More efficient ways of decomposing arrangements into topological balls have been 
proposed. In [35] the authors provide a decomposition into 0*{n'^'^~^) cells (see [iS] 
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for an improvement of this result in the case k ~ A). However, this decomposition 
does not produce a cell complex and is therefore not directly useful in computing 
the Betti numbers of the arrangement. 

8.2. Local Method. We have seen in Section [5] that in certain simple situations 
it is possible to compute the Betti numbers of an arrangement without having to 
compute a triangulation. For instance, if the arrangement has the Leray property 



(cf. Definition 5.331 Theorem 5.34 provides an efficient way of computing the Betti 
numbers o f the union. The dimension of the p-th term of the nerve complex L''(C) 
(see Eqn. (5.16 1) is this case is bounded by {p^^) = 0{nP~^^) corresponding to all 



possible {p+ l)-ary intersections amongst the n given sets. The truncated complex, 
L^_i_]^(C), can be computed by testing for non-emptiness of each of the possible 
X]i<j<£+2 (j) ~ 0{n^~^'^) at most {£ + 2)-ary intersections among the n given sets. 
The first i Betti numbers of the arrangements can then be computed from L^^j^(C) 
using algorithms from linear algebra. This technique would work, for instance, if 
one is interested in computing the Betti numbers of a union of balls in R'^. However, 
this method is no longer useful if the sets in the arrangement do not satisfy the 
Leray property. 

For non-Leray arrangements, some new ideas are needed. Before introducing 
them we first need some new notation. For the rest of this section we fix a family 

(8.1) 5 = {^i,...,5„} 

of closed and bounded semi-algebraic subsets of R'^ . For / C [n] we denote by 

(8.2) - U 

(8.3) Sj=f] S,. 

iei 

The main new idea is that in order to compute the first £ Betti numbers of a 
non-Leray arrangement S it suffices to compute triangulations, , of the sets iS'''^ 
with #J < i+2. These triangulations should have a certain compatibility property 
namely - the triangulation of Sj obtained by restricting should be a refinement 
of the triangulations of Sj obtained by restricting h"^ for all J C /. 

More formally, we define 

Definition 8.1 (Adaptive Triangulations). An ^-adaptive triangulation, h£{S), of 
5 is a collection {h^}ic[n],#i<e+2 of semi-algebraic triangulations 

(8.4) -.K' ^ S' 

having the following properties. 

(1) For each / C [n] with < £ + 2 the triangulation respects the sets 
Si,i E I. In particular, induces a triangulation of Sj, which we denote 
by hj : Kj — > Sj, where Kj is a subcomplex of . 

(2) For each ,/ C / C [n] with #/ < £ -t- 2, the triangulation hi is a refinement 
of the triangulation hj\s,. 

We now show how to obtain from a given ^-adaptive triangulation, a cell complex 
whose first £ cohomology groups are isomorphic to those of S*'"!. We will use the 
notion of homotopy colimits introduced in Section |5.9[ 
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Given an i'-adaptive triangulation, hi{S), we associate to it a cell complex, ICi(S) 
(best thought of as an infinitesimally thickened version of hocolim<^(iS)), whose 
associated topological space is homotopy equivalent to |hocolim<£(iS)|. 

Definition 8.2 (The cell complex ICi{S)). Let C denote the cell complex C(skf(Ar„i)) 



defined previously (see Definition 7.13 replacing Aq by Ar„i). Let D be a cell of 



C(sk^(Ar„i)). Then, D C |A/| for a unique simplex A/ <E Ar„i with < £ + 2 and 



(following notation introduced before in Definition 7.131 

D = Da,^,Ai n • • • n Daj^Ai n Daj, 

with /i C I2 C • • • C /p C Ip+i = I and p < £ + 1. We denote 

(8.5) IC{D) ^{Dx MM) I aeK', with hj{\(T\) C S^'}, 
and 

(8.6) JC,{S)= U IC{D). 

nec(skf(A[„j)) 

Notice that |/C£(5)| is a closed and bounded semi-algebraic set defined over 
R' = R(eo, . . . ,£f), and it contains the semi-algebraic set Ext(|hocolim<^(5)|, R'). 
Furthermore, we have 

Proposition 8.3. The semi- algebraic set 

\}Ci{S)\ 

is homotopy equivalent to 

Ext(|hocolim<<;(5)|,R'). 
Proof. From the definition of the complex JCi{S) it follows easily that 

hm|/Cf(5)| = |hocolim<^(5)|. 

It now follows (see [22, Lemma 16.17]) that |/C^(iS)| is homotopy equivalent to 
Ext(|hocolim<<?(5)|,R'). □ 



By Theorem 5.56 in order to compute the first £ Betti numbers of S*'"', it 



suffices to compute the first £ Betti numbers of |hocolim<£(5)|. Moreover, by virtue 



of Proposition 8.3 and Proposition 5.29 (homotopy invariance of the cohomology 



groups) we have that in order to compute the Betti numbers of |hocolim<£(5)| it 
suffices to compute the Betti numbers of the set |/C£(iS)|. This is the main idea 
behind the following algorithm. 

8.3. Algorithm for Computing the Betti Numbers of Arrangements. We 

can now describe our algorithm for computing the first £ Betti numbers of the set 

Algorithm 1. 

Input A family S = {Si, . . . , Sn} of closed and bounded semi-algebraic sets of R'^ 
of constant description complexity. 
Output 6o(S["l), . . . , fo^S''"!)- 
Procedure 

Step 1 Using the algorithm implicit in Theorem |4 . 5| compute an i?- adaptive triangu- 
lation, hi{S). 
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Step 2 Compute the matrices corresponding to the differentials in the co-chain com- 
plex of the the cell complex /C^(5). 

Step 3 Compute using standard algorithms from linear algebra for computing di- 
mensions of images and kernels of linear maps the dimensions of the cohomology 
groups of the complex C'{JCe{S)). 

Step 4 For < i < ^ output 

fo.CS*!"!) =dimff(C*(/Cf(5))). 

Proof of Correctness: The correctness of the algorithm is a consequence of 
Theorem |5.56| and Proposition |8.3[ □ 

Complexity Analysis: There are clearly at most ^ ( • ) ~ 0(n^+^) calls to the 

i ^ ^ 

triangulation algorithm. Each such call takes constant time under the assumption 
that the input sets have constant description complexity. Thus, the total number of 
algebraic operations (involving the coefficients of the input polynomials) is bounded 
by 0{n^~^'^). Additionally, one has to perform linear algebra on matrices of size 
bounded by 0{n^+^). □ 

9. Open Problems 

We list here some interesting open problems some of which could possibly be 
tackled in the near future. 

Computing Betti Numbers in Single Exponential Time ? Suppose S C R'" is a semi- 
algebraic set defined in terms of s polynomials, of degrees bounded by d. One of 
the most fundamental open questions in algorithmic semi-algebraic geometry, is 
whether there exists a single exponential (in k) time algorithm for computing the 
Betti numbers of S. The best we can do so far is summarized in Theorem 16.101 
which gives the existence of single exponential time algorithms for computing the 
first £ Betti numbers of S for any constant £. A big challenge is to extend these 
ideas to design an algorithm for computing all the Betti numbers of S. 

Are the Middle Betti Numbers Harder to Compute ? From the algorithm design per- 
spective it seems that computing the lowest (as well as the highest) Betti numbers 
of semi-algebraic sets, as well the Euler-Poincare characteristic of semi-algebraic 
sets, are easier than computing the "middle" Betti numbers. Is there a complexity- 
theoretic hardness result that would justify this fact? In certain mathematical 
contexts (for instance, the topology of smooth projective complex varieties) the 
middle Betti numbers contain all the information. Is there a complexity-theoretic 
analogue of this phenomenon that would justify our experience that certain Betti 
numbers are harder to compute than the others? 

More Efficient Algorithms for Computing the Number of Connected Components 
in the Quadratic Case ? For semi- algebraic sets in R'^ defined by £ quadratic in- 
equalities, there are algorithms for deciding emptiness, as well as computing sample 
points in every connected component whose complexity is bounded by k'-'^^^ [3ll45j. 
We also have an algorithm ["T for computing the Euler-Poincare characteristic of 
such sets whose complexity is k'-'^^K However, the best known algorithm for com- 
puting the number of connected components of such sets has complexity fc^*^* ' (as a 
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special case of the algorithm for computing all the Betti numbers given in Theorem 



2.6 1. This raises the question whether there exists a more efficient algorithm with 



complexity fc*^^^-* or even k^^^ ^ for counting the number of connected components 
of such sets. Roadmap type constructions used for counting connected components 
in the case of general semi-algebraic sets cannot be directly employed in this con- 
text, because such algorithms will have complexity exponential in k. 



More Efficient Algorithms for Computing the Number of Connected Components 
for General Semi-algebraic Sets ? A very interesting open question is whether the 
exponent 0{k^) in the complexity of roadmap algorithms (cf. Theorem 4.151 can be 
improved to 0{k), so that the complexity of testing connectivity becomes asymp- 
totically the same as that of testing emptiness of a semi-algebraic set (cf. Theorem 



4.121. 



Such an improvement would go a long way in making this algorithm practically 
useful. It would also be of interest for studying metric properties of semi-algebraic 
sets because of the following. Applying Crofton's formula from integral geometry 
(see for example [52]) one immediately obtains as a corollary of Theorem 4.15 (using 
the same notation as in the theorem) an upper bound of s'"' +irfO(fc ) ^j-^ ^j-^g length of 
a semi-algebraic connecting path connecting two points in any connected component 
of S (assuming that S is contained in the unit ball centered at the origin). An 
improvement in the complexity of algorithms for constructing connecting paths 
(such as the roadmap algorithm) would also improve the bound on the length of 
connecting paths. Recent results due to D'Acunto and Kurdyka [35] show that it 
is possible to construct semi-algebraic paths of length dP^''^ between two points 
of S (assuming that 5* is a connected component of a real algebraic set contained 
in the unit ball defined by polynomials of degree d). However, the semi-algebraic 
complexity of such paths cannot be bounded in terms of the parameters d and k. 
The improvement in the complexity suggested above, apart from its algorithmic 
significance, would also be an effective version of the results in [55] . 
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